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Recent  advances  in  semiconduct m technology  have  made  microprocessor 
chips  and  their  associated  interface  and  peripheral  chips  available  at 
low  cost  and  In  readily  adaptable  configurations.  1'he  availability 
ot  these  powerful  vet  small  devices  opens  up  whole  new  fields  ot 
computer  control.  1'hls  project  investigated  the  use  of  microprocessors 
and  associated  hardware  to  control  the  heating  ol  a building  in  order 
to  save  energy.  A house  was  simulated  on  the  hybrid  computer  and 
controlled  by  a microprocessor-based  digital  controller.  I’he 

microprocessor  was  a commerc i a l 1 v available  model,  one  which  could 

* 

be  obtained  at  low  cost  and  would  be  accessible  to  any  home  owner 
By  turning  the  heat  off  in  rooms  which  were  not  to  be  used,  a 
considerable  amount  of  energy  could  be  saved.  The  program  to  control 
the  room  t emperat vtres  would  be  determined  and  written  on  the  basis  ot 
the  probable  use  of  rooms.  Tests  showed  that  a significant  amount  of 
energy  was  saved  by  using  the  controller,  as  much  as  SO’.,  depending 
on  the  use-habits  ot  the  house.  The  house  simulation  can  also  be 
used  to  study  such  things  as  insulation  con  I igurat  ion , heater  sire 
and  location,  and  room  shapes  and  sizes  with  respect  to  the  amount  ot 
energy  saved  .^^Add i t ional lv , the  small  amount  of  time  used  by  the 
microprocessor  \n  controlling  the  house  suggest s the  possibility  of 

controlling  othefc  aspects  of  the  house  in  a time-sharing  modi' . The  low 

1 

cost  and  easily  configured  nature  ot  the  microprocessor-based  controller 
make  the  savings  in  energy  obtained  more  than  enough  to  justify  tlto 
installation  of  the  device  in  buildings  with  easily  adaptable  heattng 
systems.  I , JT"!  ’ T 
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PREFACE 

The  United  States  and  the  entire  world  are  approaching  a point 
where  energy  consumption  levels  must  be  drastically  reduced  or  there 
will  be  social,  economic,  and  political  changes  the  magnitude  of 
which  we  have  not  seen  since  the  great  depression  of  the  I930's. 
Economists  and  energy  analysts  have  placed  the  date  for  this  impending 
crisis,  by  careful  and  conservative  extrapolation,  at  about  the  vear 
2000.  This  is  less  than  twenty-two  years  away  from  today.  The 
massive  shortages  encountered  during  the  oil  embargo  of  1973-1974, 
the  rapidly  rising  cost  of  electrical  power,  gasoline,  oil  and  coal, 
the  recent  coal  strikes  and  their  tremendous  effects  on  the  electrical 
power  industry,  and  the  ever-diminishing  supplies  of  domestic  oil 
and  coal  all  point  toward  this  crisis.  1 feel  that  it  is  the  duty 
and  obligation  of  every  person  in  this  country  to  make  every  effort 
to  conserve  and  use  properly  the  energy  sources  we  now  have  available. 
The  President  of  the  United  States  expressed  this  need  by  declaring 
"the  moral  equivalent  of  war"  on  the  impending  crisis. 


To  fulfill  my  part  in  this  campaign  I decided  to  use  my  time 
and  effort  to  develop  a system  which  may  save  a significant  percentage 
of  the  energy  used  to  heat  homes,  buildings,  and  other  inhabited 
structures.  The  Trident  Scholar  program  gave  me  the  opportunity 
to  accomplish  this  task  while  completing  my  undergraduate  education. 
Upon  assuring  myself  that  I could  complete  the  courses  in  my  major. 
Systems  Engineering,  as  well  as  take  the  Humanities  electives  T 
consider  essential  to  any  undergraduate  curriculum  and  take  the 
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Trident  Scholar  program,  1 applied  ior  it.  My  proposal  was  to  design 
a microprocessor-based  digital  controller  which  could  be  adapted 
and  programmed  to  control  the  heating  of  a house  and  to  develop 
an  energy  saving  strategy  of  programming  to  implement  on  the 
device . 

The  great  difficulties  which  would  have  been  encountered  in 
controlling  an  actual  house  were  overcome  by  developing  a matliematical 
model  of  a house.  The  development  of  this  model  involved  techniques 
which  are  used  in  devising  the  simulation  of  any  physical  system. 

I'lie  relationship  between  the  variables  of  Interest,  the  room 
temperatures,  and  their  relationship  to  the  input,  the  outside 
temperature,  were  determined.  The  methods  used  to  find  these 
relationships  are  general  and  may  be  applied  to  any  matliematical 
model  Involving  heat  transfer  and  temperature  levels.  Some  simplifying 
assumptions  which  did  not  greatly  affect  the  accuracy  of  the  model 
were  made.  These  assumptions  simplified  the  equations  of  the 
mathemat ical  model  considerably. 


I'lie  finished  model  consisted  of  ten  simultaneous,  first  order. 


differential  equations.  This  model  was  programmed  on  the  Dartmouth 


Time  Sharing  System  using  a general  purpose,  FORTRAN  digital 
simulation  program,  OIGISIM,  developed  hy  F, . K.  Mitchell.  This 
was  done  to  get  a quick  and  easy  check  on  the  model's  validity. 
The  Initial  runs  revealed  some  errors  in  the  coefficients  causing 
the  model  to  "blow  up."  These  errors  were  corrected  and  the 


equations  were  programmed  on  the  EAI-681  Analog  computer  which. 

In  conjunction  with  a POP-IS  minicomputer  and  associated  digital-to 
analog  converters,  analog-to-digital  converters,  and  software, 
constituted  the  house  simulation.  The  general  nature  of  the 
simulation  cannot  be  stressed  enough.  These  simulation  techniques 
can  be  applied  to  any  building  to  determine  such  factors  as  the 
amount  of  energy  saved  by  changing  insulation,  adding  storm  windows, 
and  other  modifications.  These  building  modifications  would  require 
a large  amount  of  money  and  time  to  be  implemented  in  a real  house. 

In  this  simulation  they  can  be  realized  by  simply  re-computing  the 
coefficients  and  changing  potentiometer  settings.  The  simulation 
has  the  disadvantage  of  being  only  a simulation  and  therefore 
having  only  a limited  relation  to  reality.  With  care  and  good 
engineering  approximations  however,  the  results  will  closely 
simulate  the  response  of  the  real  structure.  I feel  that  the 
development  of  this  computer  simulation  was  a very  valuable  side 
product  of  my  Trident  Scholar  project  and  its  use  in  developing 
more  energy  efficient  structures  has  great  potential.  I welcome 
the  use  of  this  program  or  the  theory  used  to  develop  it  by  anyone 
who  might  be  pursuing  the  research  of  energy  efficient  house 
construction.  I feel  that  a good  deal  of  time  and  money  could 
be  saved  through  the  use  of  this  program. 

In  the  development  of  the  digital  controller,  a multitude  of 
problems  were  encountered.  I had  almost  no  background  in  micro- 
electronics or  microprocessors  and  their  construction  and  programming. 
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After  consulting  many  manuals  and  specification  sheets  and  asking 
many  questions  I was  able  to  accomplish  the  necessary  tasks.  A 
good  deal  of  the  answers  to  mv  questions  and  much  help  came  from 
my  advisors.  Professor  E.  E.  Mitchell  and  Major  Richard  Kopka . 

Bill  Lowe  and  the  hard  working  Technical  Support  Division  gave  me 
thousands  of  answers,  hints,  and  parts  as  well  as  keeping  the 
Hybrid  Computer  facilities  ready  and  usable. 

During  the  course  of  the  year  I have  learned  a tremendous  amount 
about  microprocessors,  digital  control,  simulation,  electrical 
engineering,  thermodynamics,  and  other  fields  1 would  not  have 
been  exposed  to  otherwise.  The  formal  courses  1 previously  pursued 
in  Sampled  Data  and  Digital  Control,  Electrical  Engineering,  Thermo- 
dynamics, and  Computer  Simulation  were  all  of  great  use  and  the 
knowledge  I had  retained  from  these  courses  was  reinforced  and 
supplemented  by  its  application  in  the  project.  The  Trident  Scholar 
program  was  the  most  rewarding  means  of  completing  my  four  years 
at  the  Naval  Academy,  both  in  terms  of  doing  my  part  in  the  war 
on  excessive  energy  consumption  and  in  preparing  myself  for  service 
as  an  Officer  in  the  Navy. 
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CHAPTER  1 
IN 1R0DUCT10N 

l'he  basic  purpose  of  tills  project  was  tlic  digital  control  ol  an 
environmental  svstem  such  as  a house  or  building  where  a specific 
environment  must  be  maintained  by  the  addition  or  subtraction  of 
heat  to  individual  rooms,  groups  of  rooms,  or  the  whole  structure. 

Due  to  the  possibly  superfluous  nature  of  air  cond i t ion ing (cool ing) 
and  the  possibility  that  its  use  will  be  curtailed  in  any  reasonable 
emergency  energy  program,  the  controller  has  been  designed  to  control 
only  the  heating  of  a structure.  The  use  of  this  controller  to 
control  both  the  heating  and  cooling  of  a structure  is  entirely 
possible,  with  some  additional  programming  it  could  be  used  satis- 
tactorily  to  control  the  cooling  of  a structure.  The  control 
concepts  to  be  considered  are  those  involving  a microprocessor 
controlled  sampler  and  decision  making  device.  The  use  of  an  actual 
structure  to  test  and  evaluate  the  controller  was  not  feasible  due 
to  the  typical  failure  of  outside  temperatures  to  follow  test  parameters 
and  the  large  cost  of  rigging  a building  for  alternate  control  bv 
both  digital  controller  and  conventional  thermostats.  Therefore 
an  ana log/hvbrid  computer  simulation  of  a building  was  used  to 
test  and  develop  the  digital  controller.  The  ease  of  varying 
input  parameters,  access  to  measurable  outputs,  and  speed  of 
test  runs  are  some  of  the  advantages  of  using  an  analog/hybrid 
computer  simulation  as  opposed  to  the  use  of  actual  physical 
structures.  Other  benefits  include  the  ability  to  actually  model 
the  entire  control  system,  both  the  digital  controller  and  the 
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convent  Iona  1 thermostat  control  on  the  digital  portion  ot  the  hvhrtd 
Computer.  Thus  the  Ideal  performance  of  the  digital  controller  can 
he  measured  before  It  is  built,  giving  some  indication  ot  its 
appl  tcab  1 1 1 tv  . The  ability  to  develop  a suhopt  Imum  programmi  ng 
routine  for  the  microprocessor  using  the  hybrid  computet  made  the 
implementat ion  ot  the  digital  controller  much  easier  than  it  all 
programming  had  been  developed  on  the  microprocessor  using  the 
appropriate  machine  language . 

The  protect  was  carried  out  in  a series  of  stages,  each  building 
on  the  last  and  concluding  with  tests  of  the  amount  of  energy  that 
could  be  actually  saved  by  vising  the  controller.  The  first  stage 
was  the  development  of  a hybrid  computer  simulation  of  a house  and 
surrounding  environment  with  respect  to  heat  flow  and  temperature 
levels.  Certain  assumptions  were  made  tv’  simplify  this  simulation 
and  these  will  be  discussed  in  detail  later.  The  second  stage  was 
to  develop  both  conventional  and  digital  control  algorithms  on  the 
digital  portion  ot  the  hybrid  computer.  The  digital  control  algorithm 
was  simplified  for  ease  of  implementation  ot  the  programs  on  the 
microprocessor , due  to  the  difficult  nature  of  machine  language 
programming.  The  third  stage  of  the  project  was  to  build  the 
microprocessor-based  computer  and  an  interlace  device  to  connect 
the  computer  to  the  simulation.  The  microprocessor  used  was  the 
Intel  8080,  as  configured  in  the  Heath  H-8  digital  computer.  The 
Heath  H-8  was  chosen  because  ot  its  easily  expandable  memory , the 
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interface  capabi  1 it v , and  the  fact  that  this  unit  was  easily  obtained 
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and  assembled  bv  the  user.  The  implementation  of  the  digital  controller 
would  be  within  the  grasp  of  the  common  user  rather  than  limited  to 
the  skilled  technician.  The  microcomputer  was  then  linked  up  with 
the  simulation  and  the  control  algorithms  were  programmed  on  the 
machine.  The  final  stage  of  the  protect  was  to  run  many  tests  to 
evaluate  the  performance  of  the  digital  controller  with  respect  to 
the  amount  of  energy  saved  when  compared  to  the  energy  used  bv  conventional 
thermostat  control  of  the  heating  of  the  structure.  The  performance 
of  the  device  will  give  some  indication  of  the  feasibility  of  implementing 
it  on  existing  houses  and  of  building  the  control  device  into  a new 
h ome . 
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A MU  HI  MU  U AL  MODI. I . 


The  simplification  of  the  control  problem  brought  about  by 
ci  -itrol 1 ing  a simulation  rather  than  an  actual  structure  led  to 
the  development  of  a mathematical  model  of  a tvpical  house.  A 
five  room,  one  story  house  with  closets,  two  bedrooms,  kitchen 
and  living  room  was  used  as  the  modeled  structure.  A sloping 
roof  and  attic  were  also  included  in  the  house  model (see  Appendix  A1 


The  materials  used  in  the  construction  of  walls,  doors,  windows, 
ceilings,  and  floors  were  all  chosen  as  typical  of  many  homes.  In 
addition,  the  fact  that  this  is  a mathematical  model  made  the 
rapid  substitution  of  different  materials  possible,  by  substituting 
the  new  heat  transfer  coefficient  into  the  governing  equation.  This 
enabled  the  effects  of  changing  insulation  to  be  quickly  determined, 
and  indications  of  the  feasibility  of  adding  or  changing  insulation 
in  a house  to  be  quickly  obtained.  The  simulation  was  also  easily 
changed  to  represent  different  types  of  houses,  the  onlv  limiting 
factor  being  the  number  of  non-cont iguous  rooms  in  the  house.  The 
analog  computer  used  bad  a limited  number  of  integrators  and  therefore 
could  only  be  programmed  to  solve  an  equally  limited  number  of 
differential  equations.  The  possibility  of  using  larger  computers 
or  slaving  several  of  the  computers  together  to  achieve  simulation 
of  larger  houses  or  buildings  is  obvious  and  bears  further  invest i ga t i o« 


At  this  point  the  basic  governing  equations  were  used  to  define 


describing  the  temperature  of  a room  become  very  large  and  increasingly 
difficult  to  program.  To  avoid  this  needless  comp  1 icat ion , and  to 
get  a clearer  picture  of  the  effects  of  the  control  on  the  variables 
of  interest,  certain  simplifying  assumptions  were  made.  The  first 
assumption  made  is  that  of  constant  air  density.  The  value  of  the 
air  density  inside  the  house  was  assumed  to  be  a constant  value  of 
.0764  pounds  per  cubic  foot.  The  second  assumption  is  that  all 
doors  and  windows  remain  closed.  This  would  be  the  case  in  very 
cold  weather  when  doors  would  be  opened  for  a minimum  of  time  and 
windows  would  be  kept  closed.  The  third  assumption  is  that  of 
perfect  ventilation  and  convection.  The  air  in  each  room  is  assumed 
to  be  at  a constant  temperature  throughout  that  room.  The  temperature 
of  a room  is  the  same  regardless  of  the  location  within  the  room. 

The  fourth  assumption  is  that  the  infiltration  effect  is  considered 
negligible,  that  no  wind  infiltrates  through  the  walls  or  the 
windows . 

1 

All  these  assumptions  are  not  really  a factor  in  determining 
tlie  results  which  are  desired  of  the  simulation.  The  goal  is  to 
control  the  heating  of  rooms  in  such  a way  that  heat  is  conserved 
when  compared  with  conventional  methods  of  controlling  the  heating 
of  houses.  When  assumptions  are  made  which  are  typical  of  well 


i t i 


built,  heat  conservative  houses,  the  task  of  saving  energy  becomes  a 
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more  difficult  procedure  and  the  results  will  he  more  significant 
than  the  saving  of  energy  in  a low  efficiency  house  would  he.  Bv 
modeling  a house  which  is  already  energy  efficient,  the  task  of 
controlling  the  heating  of  this  house  in  a manner  which  saves  energy 
becomes  a much  more  significant  and  important  task. 

The  differential  equations  which  relate  the  temperature  of  a 
room  to  the  temperature  of  surrounding  spaces  is  derived  f rom  two 
basic  laws  of  heat  flow.  The  lirst  relates  the  heat  flow  across  a 
heat  conducting  boundary,  dQ/dt,  to  the  coefficients  of  heat  transfer 
in  that  boundary,  k,  the  temperature  differential  across  that  boundary, 
dT,  the  thickness  of  the  boundary,  d.X,  and  the  area  of  the  boundary, 

A (see  Appendix  A). 

dQ/dt  = dT/dX  * k * A 

The  change  of  temperature  of  a thermal  mass,  in  this  case  a room,  is 
related  to  the  thermal  capacity  of  the  material,  Cp*W,  and  the  hear 
flow  into  the  material.  W is  the  weight  of  the  air  in  the  room  and 
Cp  is  the  thermal  constant  of  the  air. 

dT  = l/(Cp  * W)  * dQ 

These  equations  relate  the  temperature  inside  a room  to 
the  temperature  of  the  surrounding  spaces.  The  method  used  to  get 
the  equations  for  each  room  in  the  house  to  be  simulated  is  not  trivial 
so  the  process  will  he  considered  in  detail  for  one  of  the  rooms. 

The  heat  transmission  coefficients  which  combine  the  thermal  conductivity 
of  the  materials,  the  thickness  of  the  walls  and  the  effects  of  such 
things  as  the  film  effect  and  other  factors,  are  obtained  from  a table 
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of  typical  values  for  walls,  doors,  windows,  floors,  and  ceilings. 

See  Appendix  A for  a table  of  the  values  used  In  this  simulation. 

These  lumped  coefficients  allow  the  equation  formulation  process 
to  proceed  with  much  more  ease  than  if  the  walls  and  other  boundaries 
were  taken  material  by  material  and  the  exact  thermal  transmission 
coefficient  were  obtained  for  each  boundary.  The  living  room  (room  Al 
is  measured  and  areas  of  the  interface  with  spaces  adjacent  are  found. 

Meat  flow  is  the  product  of  the  lumped  transmission  coefficient,  u,  the 
area  of  the  wall,  and  the  temperature  differential  across  the  wall. 

The  interface  with  the  outside  consists  of  two  walls,  three  windows, 
and  a door.  The  respective  areas  are: 
walls  = 195  square  feet 
windows  = 18  square  feet 
door  =•  19  square  feet 

When  multiplied  bv  the  appropriate  u factors  the  total  coefficient 
of  the  outside  temperature  becomes: 

2 9 9 

Heat  flow  » (Tout  - Ta)(195  ft  * .24  BTU/hr-ft  -dog . F + 19  ft 

* .48  BTU/hr- ft^-deg . F + 18  ft2  * 1.243  BTU/hr- ft  -deg . F) 

Heat  flow  - (Tout  - Ta)(78.293  BTU/hr-deg.  F> 

The  heat  flow  through  the  interface  between  room  A and  room  E is  computed. 

2 

interior  walls  » 10.33  ft 

2 

Interior  door  « 19  ft 

Heat  flow  - (Te  - Tal (10.33  ft2  * .33  BTU/hr- ft2-deg . F + 

19  ft2  * .64  BTU/hr- ft2 -deg.  F) 

Heat  flow  - (Te  - Tal (15.57  BTU/hr-deg.  Fl 


The  heat  flow  through  the  Interface  between  room  A and  room  F Is 

2 

Interior  walls  « 63.6b  It 

9 

interior  door  ■ 19  ft 

Heat  flow  - (Tf  - Ta) (63.66  ft2  * .33  BTU/hr-ft2-deg . F + 

19  ft2  * .64  BTU/hr-ft2-deg . F) 

Heat  flow  = (Tf  - Ta) (33.17  BTU/hr-deg.  F) 

The  heat  flow  through  the  interface  between  room  A and  room  H is 

interior  wall  = 29  ft2 

2 

interior  door  » 19  it 

Heat  flow  = (Th  - Ta)(29  ft2  * .33  BTU/hr- ft2-deg . F + 19  f 
.64  BTU/hr- ft2 -deg.  F) 

Heat  flow  = (Th  - Ta)(21.73  BTU/hr-deg.  F) 

The  heat  flow  through  the  interface  between  room  A and  room  I!  is 
interior  wall  = 72  ft2 

Heat  flow  = (Tb  - Ta)(72  ft2  * .33  BTU/hr- ft2-deg . F) 

Heat  flow  = (Tb  - l'a)(23.76  BTU/hr-deg.  F) 

The  heat  flow  through  the  interface  between  the  ceiling  and  room 
computed . 

2 

area  of  ceiling  = 176  ft 

Heat  flow  = (Tci  - Ta)(176  ft2  * .12  BTU/hr- ft2 -deg . F) 

Heat  flow  = (Tci  - Ta) (21.12  BTU/hr-deg.  F) 


compu  t ed 


computed 


2 * 


computed 


The  heat  flow  through  the  interface  between  the  floor  and  room  A 


IS 


is  computed. 

2 

Area  ol  floor  “ 176  ft" 

Heat  flow  - (Til  - Ta)(176  ft'  * .!>8  BTU/hr- i t" -deg . 1") 

Heat  flow  = (Tfl  - la)  (102. 08  BTU/hr-deg.  F) 

rhe  total  heat  flow  into  room  A is  the  sum  ol  all  these 
individual  heat  flows,  plus  any  heat  source  or  sink(heater  or  cooler). 
dQ/dt  = 78 . 2 94 (Tout  - la)  + 15.57(Te  - Ta)  + 33 . 1 7 (If  - Ta)  + 

2 1 . 73  (111  - Ta)  + 2 3.76  (lb  - la)  + 102.08(111  - Ta)  + 21.12(Tci  - Ta) 
+ Qexternal (heater) 

When  all  terms  are  multiplied  out  and  the  common  terms  are  collected, 
the  equation  becomes: 

dQ/dt  *=  -295.724  * Ta  + 78.294  * Tout  + 1S.S7  * Te  + 33.17  * Tf 
+ 21.73  * lh  + 23.76  * Tb  + 102.08  * Tfl  + 21.12  * Tci  + Qext . 

This  equation  gives  the  heat:  flow  into  room  A as  a function  of  the 
temperatures  of  adjacent  spaces  and  the  applied  heater.  The  units 
are  BiTI/hr.  A quick  check  on  the  validity  of  the  equation  is  to 
let  all  temperatures  except  that  in  room  A be  set  equal  to  that  in 
room  A.  When  this  is  done,  the  heat  flow  equation  becomes: 
dQ/dt  = Qext . 
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Ibis  equation  is  useful  but  what  is  required  is  a relation  between 
the  temperature  in  room  A and  the  temperatures  of  adjacent  spaces. 

To  find  this  relation  the  thermal  mass  of  room  A must  he  computed. 

The  thermal  mass  in  the  simulation  is  assumed  to  be  only  the  air 
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lb 

within  the  room.  This  mav  not  be  an  entirely  accurate  assumption, 
but  the  only  thins  affected  by  the  changing  of  the  thermal  mass  is 
the  time  constant  of  the  room.  The  steady  state  response  is  not 
changed  and  the  total  amount  of  heat  is  not  changed.  We  must  now 
find  the  weight  of  the  air  in  the  room  considered,  room  A.  This 
is  done  by  finding  the  volume  of  the  room  and  multiplying  It  bv 
the  density  of  air. 

Volume  of  room  A = 17b  ft"  * 8 ft  « 1408  ft^ 

3 -» 

W = weight  of  the  air  = .074b  lhs/ft  * 1408  ft  = 107.5712  lbs. 

1/CpW  ■ (1/.24  BTU/lbs-deg  Fl (1/107.5712  lbs!  = .03873  deg  F/BTU 
This  factor  relates  the  heat  flow  into  room  A to  the  temperature  in 
room  A.  When  the  equation  for  the  heat  flow  into  room  A is  multiplied 
by  this  factor,  the  result  is: 

dT/dt  = -11.454  * Ta  + 3.033  * Tout  + .6031  * Te  + 1.285  Tf  + 

.8417  * Th  + .9203  * Tb  + 3.954  * Tfl  + .8181  * Tci  + .03873  * Qext . 
(deg.  F/hrl 

This  equation  relates  the  temperature  in  room  A to  the  temperature 
in  adjacent  spaces  but  the  coefficients  are  not  in  the  best  form  for 
programming  on  an  analog  computer  so  the  change  in  temperature  per 
minute  is  found  by  dividing  the  eqvtat  ion  bv  60. 

dT/dt  = -.1909  * Ta  + .0505  * Tout  + .0101  + Te  + .0124  Tf  + 

.014  * Th  + .0153  * Tb  + .0659  * Tfl  + .00646  * Qext  + .0136  * Tci. 

(deg  K/mtnl 

The  same  procedure  is  followed  for  all  ten  rooms  to  arrive  at  the 
complete  model.  There  are  ten  compartment  temperatures,  five  heaters, 
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and  the  outside  temperature  as  variables  in  the  equations.  The  outside 
temperature  and  the  heaters  are  considered  inputs  and  the  temperatures 
are  found  by  solving  simultaneously  the  ten  first  order  differential 
equations  which  comprise  the  model.  See  appendix  A for  a complete 
list  of  the  equations  and  the  procedure  used  for  each  room  to  arrive 
at  these  equations. 

This  model  represents  a frame  house  upon  which  certain  restrictions 
have  been  placed.  While  this  may  not  give  a totally  accurate  picture 
of  the  characteristics  of  the  house,  the  factors  of  interest  are 
simulated  with  a high  degree  of  correlation  with  the  actual  house, 
assuming  the  restrictions  are  true  of  the  real  house  as  well.  The 
simulation  which  includes  every  variable  and  relation  that  could 
possiblv  affect  the  parameters  of  interest  mav  very  well  be  impossible 
to  program.  If  it  is  possible,  the  program  will  be  quite  large  and  the 
result  ing  view  of  the  parameters  of  interest  will  be  fuzzy  and 
probably  inaccurate  due  to  the  accumulation  ol  the  small  amounts 
of  machine  error,  static  charge,  and  other  stochastic  factors  which 
are  a part  of  every  machine  realized  simulation.  A simulation  which 
neglects  the  minor  relations  and  considers  the  ma ior  relations  gives  a 
clear,  more  easily  realized  view  of  the  parameters  in  question.  Kven 
though  the  results  mav  be  slightly  different  from  the  real  svstem,  the 
use  of  computer  simulation  to  accurately  and  quicklv  predict  the 
performance  of  the  actual  system  is  a proven  fact  . l'his  model  of  a 
house  can  he  used  to  predict  the  behavior  of  an  actual  house  which 
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is  exposed  to  varying  outside  temperatures. 
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l'he  simulation  now  consists  ot  ton  simultaneous  differential 
equations  with  ton  temperature  variables,  five  heater  inputs,  an 
outside  temperature  input,  and  the  floor  temperature  input.  To  solve 
these  equations  bv  hand  would  he  very  difficult  unless  the  inputs 
were  all  assumed  to  be  constants  and  the  temperatures  were  also 
assumed  to  be  constants.  Then  the  derivatives  would  all  he  equal 
to  zero  and  the  solution  becomes  a matter  of  solving  ten  simultaneous 
algebraic  equations.  l'he  answer  to  the  problem  of  solving  these 
equations  in  such  a manner  as  to  make  the  simulation  useful  is  to 
mechanize  the  solution.  That  is,  use  the  computer  to  solve  the 
equations  in  a high  speed  mode,  giving  a nearlv  continuous  solution. 

1 he  entire  model  was  programmed  on  the  Dartmouth  Time  Sharing  System 
using  a general  purpose  digital  simulation  program,  DIGISIM,  developed 
bv  K.  K.  Mitchell.  This  solution  gave  some  insight  into  the  care 
originally  exercised  in  computing  the  coefficients  of  the  equations, 
as  the  simulation  "blew  up"  when  run.  Alter  recomputing  the  equations 
and  correcting  the  errors,  the  house  followed  the  outside  temperature 
as  it  should,  with  a time  lag  dependent  on  the  room.  l'he  controller 
was  also  Included  with  this  program  and  appeared  to  be  functioning  well. 
When  the  heat  used  over  a period  of  time  was  measured  with  the  rooms 
maintained  at  some  temperature  above  the  outside  temperature,  a 
change  was  noted  when  the  delta  time  input  to  the  program  was  changed. 
Attempts  to  make  the  delta  time  Increment  as  small  as  was  necessary 
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to  assure  accurate  solution  bv  the  Kunge-Kutta  method  resulted  in  verv 
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long  run  times  for  the  simulation  without  assuring  a good  solution. 
The  heat  used  over  a period  of  time  was  seen  to  he  a function  of  the 
time  increment  used  to  solve  the  Runge-Kutta  integration  algorithm. 

A listing  of  this  program  is  included  in  Appendix  H. 
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The  simulation  was  moved  to  the  analog  eompul  ur  where  the 
dynamics  were  programmed  on  the  KAl-hHl  analog  computer  and  the  controller 
and  other  programs  were  programmed  on  the  l’DP-15  digital  computer.  See 
Appendix  i-  for  the  analog  l low  charts  and  listings  of  the  digital 
programs  used  in  the  simulation. 
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CHAi'Ti  K 4 
fhL  INTERFACE 


To  control  the  analog,  simulation  of  the  house  with  the 
H-8  microcomputer,  an  interface  between  the  H-8  anti  the  KAI-681 
was  necessary.  There  were  two  primary  reasons  for  this  interface. 

The  first  was  brought  about  by  the  fact  that  the  H-8  performs 
operations  on  8-bit  binary  numbers,  while  the  EAl-881  uses  analog 
voltage  levels  to  represent  variables  in  its  solution  of  equations; 
thev  represent  room  temperatures  in  this  case.  To  convert  these 
various  voltages  into  information  which  could  be  processed  by  the 
H-8  in  order  to  control  the  simulation,  an  interface  was  needed.  The 
second  reason  for  an  interface  was  the  incompatibility  of  the  H-8 
control  outputs  and  the  voltages  needed  to  control  electronic 
switches  on  the  EAl-681.  Also  the  possibility  of  introducing  a 
voltage  level  which  is  not  of  the  type  expected  by  the  digital 
system,  such  as  a negative  voltage  level  existed.  So  a buffer 
svstem  between  the  digital  output  of  the  H-8  and  the  digital 
portion  of  the  F.Al-681  was  needed. 

An  8-blt  analog-to-dtgital  converter  was  used  to  convert  the 
analog  voltage  levels  into  binary,  8-bit  numbers  for  the  H-8  to 
process.  There  were  several  rooms  to  be  monitored  so  there  was  a 
need  for  either  several  analog-to-digltal  converters  or  multiplexed 
input  to  a single  analog-to-dlgital  converter.  Since  the  H-8  could 
onlv  read  one  parallel  I/O  port  at  a time,  receiving  one  8-bit 


number,  then  another,  a multiplexer  would  have  been  necessary  to 
select  which  analog-to-digital  converter  to  read.  To  save  space, 
time  and  materials,  a multiplexer  was  used  to  select  one  signal 
to  go  into  a single  analog-to-digital  converter  which  was  read 
bv  a parallel  I/O  port  on  the  H-8  when  the  "data  valid"  signal 
went  low.  The  multiplexer  was  controlled  by  the  output  from  one 
of  the  parallel  I/O  ports  and  the  data  from  the  analog-to-digital 
converter  was  read  into  the  input  side  of  the  same  port.  To 
protect  the  analog-to-digital  converter  from  the  possible  spurious 
voltage  levels  that  might  occur  in  the  KA1-681,  operational  amplifier 
voltage  followers  were  used  as  analog  buffers  on  the  inputs  to  the 
multiplexer  and  on  the  input  to  the  analog-to-digital  converter. 


The  basic  component  for  the  analog-to-digital  interface  was 
an  ADC-EK8R  monolithic  analog-to-digi tal  converter  chip  manufactured 
bv  Patel  Systems  Inc. (see  specification  sheet  11.  The  input  for 
tills  chip  came  from  a 00401 1 A analog  mult  ip  lexer (see  specification 
sheet  21  through  a |J74l  operational  amplifier  voltage  follower 
(see  specification  sheet  41.  The  multiplexer  was  controlled  bv 
TTL  signals  from  the  H-8  which  were  pulled  up  to  the  logic  level 
needed  to  operate  the  CH4011A.  These  signals  were  pulled  up  to 
the  necessary  level  bv  putting  them  through  a N7407  non-inverting 
buf fer/drlver (see  specification  sheet  11.  This  buffer  has  open 
collector  outputs,  allowing  the  signals  to  be  pulled  up  with  a 
voltage  source  and  a resistor  for  each  signal.  The  output  of  the 
analog-to-digital  converter  chip  was  buffered  with  a CP4010AE 
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hex  bu f fer/ converter ( see  specification  sheet  3).  These  signals 
then  go  into  the  parallel  I/O  port  on  the  H-8  where  they  are  read 
by  the  CPU.  Control  for  the  ADC-EK8B  analog-to-digital  converter 
was  accomplished  in  an  indirect  manner.  The  signal  which  told  the 
chip  to  "start  convert"  was  connected  to  the  "send  data"  output  of 
a parallel  I/O  port  on  the  H-8.  This  signal  was  high  all  the  time 
except  for  a pulse  of  about  five  microseconds  duration,  occuring 
whenever  the  port  is  read  by  the  CPU.  This  is  because  of  the  inverse 
logic  used  by  Heath  peripherals.  Therefore,  the  analog-to-digital 
converter  was  essentially  in  a continuous  mode(see  signal  relationship 
chart).  The  "data  valid"  signal  from  the  analog-to-digital  converter 
was  sent  via  the  CD4010AE  buffer  to  the  "data  sent"  input  on  the 
parallel  I/O  port.  The  hardware  of  the  port  insured  that  the  port 
would  not  be  read  unless  the  signal  present  at  the  "data  sent"  input 
was  making  a transition  from  high  to  low.  Since  the  chip  was  converting 
continuously,  the  changing  voltage  levels  would  be  read  at  the  value 
they  were  when  the  I/O  port  was  read,  not  the  value  they  were  at  when 
the  "start  convert"  signal  was  sent  to  the  converter  chip.  This 
helped  eliminate  some  of  the  time  lag  involved  in  controlling  the 
s imulat ion . 


The  H-8  was  programmed  to  use  this  interface  in  controlling 
the  simulation.  The  control  outputs  were  sent  from  the  parallel 
I/O  port  where  they  were  latched  between  output  commands,  into  an 
N7407  buffer.  On  the  output  side  of  the  buffer  the  signals  were 
pulled  up  to  the  TTL  logic  level  used  in  the  F.A1 -681  analog  computer. 
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positive  five  volts.  The  buffer  also  assured  that  the  output  signals 
would  have  the  necessary  fan-out  strength.  These  signals  were  used 
to  turn  electronic  switches  on  and  off,  thereby  simulating  the  turning 
on  and  off  of  heaters. 

The  H-8  computer  was  able  to  control  the  simulation  by  taking 
a set  of  analog  voltages  and  outputlng  a set  of  control  signals, 
with  the  interface  making  the  signals  compatible  both  to  the  H-8 
and  to  the  EAI-681.  Now  the  H-8  can  be  programmed  to  control  the 
simulation  without  regard  for  the  basic  incompatibility  of  the  signals. 
Ideallv  the  interface  should  be  transparent  to  the  H-8  both  in  terms 
of  input  and  output.  This  quality  is  achieved  almost  perfectly  on 
the  output  of  die  H-8,  but  the  inherent  time  lag  involved  with 
converting  analog  voltages  into  binary  numbers  make  transparency 
of  the  input  to  the  H-8  impossible.  However  with  the  ADC-EK8B 
converter  in  the  continuous  mode,  the  time  lag  is  minimized. 
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ADC-KK8B  Signal  Relationship 
Chart 
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The  parallel  port  on  the  11-8  is  set  up  to  read  the  port  when 
the  "data  sent"  input  makes  a transition  from  high  to  low.  Thus  the 
ADC  will  be  read  when  the  data  valid  signal  is  making  the  transition 
from  high  to  low.  The  start  convert  signal  comes  from  the  "send 
data"  output  on  the  11-8  parallel  port.  This  signal  is  high  except 

for  a very  short  low  pulse  when  the  port  is  rend  due  to  the  inverse 
logic  of  Heath  peripherals.  This  signal  is  seen  ns  an  effectively 

continuous  high  signal  by  the  APC-KK8B,  putting  It  in  the  free 
running  mode. 
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CONTROL 


To  have  the  11-8  computer  control  the  house  simulation,  it  was 
necessary  to  write  o program  which  could  read  the  time  and  then  the 
temperature  of  each  heated  room,  making  a decision  to  turn  the  heater 
in  each  heated  room  either  on  or  off  according  to  the  reference  for 
that  room.  The  program  to  accomplish  this  task  will  he  supplied  with 
reference  values  which  are  associated  with  a particuliar  time  period. 
The  user  will  eventually  program  both  the  parameters  of  the  time 
periods  and  the  temperature  references  for  each  room  associated  with 
that  time  period.  For  the  sake  of  simplicity  the  time  periods  were 
chosen  and  built  into  this  program  and  the  temperature  references 
were  left  as  inputs  in  the  form  of  octal  numbers,  each  of  which 
corresponds  to  a particular  voltage  level  or  temperature  level. 

It  is  necessary  for  the  program  to  select  the  multiplexer  command 
to  put  the  time  signal  into  the  analog-to-digita l converter  of 
the  interface.  The  value  which  results  is  used  to  select  which 
time  frame  and  therefore  which  set  of  temperature  reference  levels 
are  used.  The  program  then  reads  each  room  temperature  in  a similar 
manner,  compares  it  to  its  reference  level  and  adds  the  control  word 
for  that  room  to  the  command  word  or  goes  on  to  read  the  temperature 
of  the  next  room.  Since  the  control  words  are  all  in  increasing 
powers  of  two,  in  binary  numbers,  they  can  simply  be  "ORed"  together 
to  give  a total  command  word  which  will  affect  all  heated  rooms.  This 
command  word  is  then  output  to  the  control  lines  through  a parallel 
input/output  port,  turning  the  heaters  either  on  or  off.  The  program 
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then  cycles  through  again,  starting  with  the  read  time  operation.  In 
real  time  the  program  would  have  a period  of  rest  at  this  point,  but 
the  speed  of  the  simulation  made  continuous  operation  necessary.  The 
program  flow  chart,  mnemonic  code  listing,  octal  operation  code  listing, 
mnemonic  code  definitions,  and  voltage  level  to  octal  conversion  chart 
for  the  interface  ana log-to-digital  converter  are  listed  in  Appendix  E. 
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CHAPTER  0 
RESULTS 

The  house  simulation  gave  verv  satisfactory  results.  Its  behavior 
compared  favorably  with  the  performance  of  an  actual  house.  These 
conclusions  about  the  simulation  are  subjective,  as  the  house 
which  was  simulated  does  not  exist  and  any  comparison  on  a point  bv 
point  basis  was  impossible.  The  simulation  Itself  was  used  to 
develop  an  energy  efficient  programming  strategy.  In  the  project 
proposal  it  was  assumed  that  a median  temperature  setting  existed 
between  the  outside  temperature  and  the  comfortable  inside  temperature, 
which  when  vised  for  the  reference  temperature  for  an  unused  room  would 
vietd  a minimum  heat  use.  This  median  temperature  was  sought  bv 
varying  the  reference  temperature  of  one  room  while  holding  all 
others  at  an  arbitrary  temperature  with  heaters.  The  outside 
temperature  was  set  at  20  degrees  F and  all  inside  rooms  except  the 
ones  being  tested,  either  one  at  a time  or  in  combinations  of  two 
at  a time,  held  at  75  degrees  K.  The  rooms  being  tested  were 
maintained  at  various  temperatures  ranging  from  the  outside  temperature, 
20  degrees  F,  to  the  inside  temperature,  75  degrees  F.  The  simulation 
was  then  run  for  eight  hours (equivalent  real  time!  and  the  total  heat 
vised  was  measured.  The  results  of  these  test  runs  are  displayed 
graphically  in  Appendix  F.  These  tests  gave  an  indication  of  whether 
a temperature  above  the  outside  temperature  would  he  a more  heat 
efficient  setting  for  the  unused  rooms.  As  the  graphs  Indicate,  the 
heat  use  is  at  a minimum  when  the  temperature  setting  is  below  the 
steady  state  unheated  temperature  for  the  room.  The  steady  state 
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uuheated  temperature  lor  a room  is  that  temperature  which  the  room 
settles  to  when  uuheated.  The d l f ferenoe  between  this  temperature 
and  the  outside  temperature  is  caused  bv  the  heat  tlow  t rom  adjacent 
rooms.  In  all  cases  the  minimum  heat  use  occurred  when  the  unused 
room  was  uuheated.  The  temperature  that  the  room  settled  to  when 
unheated  is  a function  ot  the  room,  the  temperature  of  the  adjacent 
rooms,  and  the  outside  rooms.  These  results  are  in  accordance  with 
a theoretical  analysis  of  the  house. 
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A simplified  model  ot  a house  with  four  symmetric,  cubic  rooms, 
no  attic  and  all  walls  made  of  like  materials  is  vised  for  the  theoretical 
analysis  ot  the  house.  A simpler  structure  is  used  for  this  analysis 
because  of  the  size  and  number  ot  equations  involved  in  the  simulation 
used  in  this  project . 
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Assume  that  rooms  it,  t’,  and  P are  heated,  maintained  at  a temperature 
Ti,  greater  than  Tout.  Room  A is  uuheated  and  the  system  is  at  steady 
state,  so  that  all  derivative  terms  are  equal  to  tero. 


Ta  * (K2  * Tl  + K3  * Ti  + KS  * Tout)/KA 


The  temperature  in  room  A when  unheated  is  shown  to  he  a function 
of  the  temperatures  in  the  adjacent  rooms,  the  outside  temperature,  and 
the  construction  of  the  room.  The  heat  use  relations  between  the  rooms 
can  be  determined  in  a similar  manner.  For  the  temperature  change  to 
be  equal  to  zero,  the  net  heat  flow  into  the  room  has  to  be  equal  to 
zero.  This  means  that  the  heater  output  is  just  equal  to  the  heat 
loss  of  the  room,  either  because  of  the  size  of  the  heater  or  more 
commonly  the  heater  is  turned  off  and  on  in  such  a manner  that  the 
total  heat  over  a period  of  time  is  equal  to  the  total  heat  loss 
of  the  room.  To  accomplish  this  a controller  is  used  to  maintain 
the  temperature  in  the  room  at  a constant  level,  thereby  balancing  the 
net  heat  flow  at  zero  over  some  time  period.  To  illustrate  this 
the  simplified  model  is  again  used.  The  heat  flow  equations  for 
the  rooms  are: 


dQfAV'dt  = -KA'*Ta  + K2'*Tb  + K3’*Tc  + KV*  Tout  + Heater  A * 0 
dQ(BWdt  = -KR'*Tb  + Kl'*Ta  + K4'*Td  + K^'*  Tout  + heater  R = 0 
dQ(C)/dt  - -KC'*Tc  + Kl'*Ta  + K4'*Td  + Kr>’*  Tout  + Heater  C = 0 
dQ(Pl/dt  » -Kl)’*Td  + K2  ’*Tb  + K3'*Tc  F KV*  Tout  F Heater  H - 0 


When  these  equations  are  summed  and  the  individual  heaters  are  combined 
into  a single  heater  the  result  is: 


Heater  = (2*  Kl  ' - KA">  * Ta  -F  (2*  K2'  - KR'l  * Tb  -F  (2*  K3'  - KO ' 3 


Tc  + (2*K4 ' - KB ' 3 * Td  + 4 * K5'  * Tout 


32 

saving  of  energv.  The  onlv  consideration  Is  the  room  Itself,  whether 
It  can  be  allowed  to  go  to  the  low  temperatures  that  mav  result  when 
the  outside  temperature  is  very  low.  If  this  Is  a factor  the  room 
may  have  to  be  programmed  to  stav  above  some  low  temperature,  not  the 
most  efficient  setting  but  necessary  for  the  room.  It  has  been 
stressed  that  this  analvsis  is  valid  only  for  the  steadv  state  case. 

If  the  room  temperatures  are  changing,  the  derivative  terms  do  not 
go  to  zero  and  the  solution  becomes  a great  deal  more  complicated.  The 
solution  of  this  problem  is  one  of  the  reasons  that  the  simulation 
was  implemented  on  the  hybrid  computer. 

Tests  were  run  for  a variety  of  different  conditions  where  a 
steadv  state  solution  would  not  be  appropriate.  The  first  series 
of  tests  involving  the  setting  of  the  reference  temperature  to  a 
constant  level  and  maintaining  it  there  for  a period  of  eight  hours 
were  duplicated  with  some  changes.  The  unused  room  or  rooms  were 
initially  heated  to  the  same  level  as  the  used  rooms  and  then  turned 
to  a lower  reference  temperature  for  a period  of  time  less  than  the 
total  run  time.  The  unused  room  was  again  turned  back  up  to  the 
reference  temperature  of  the  used  rooms  before  the  test  run  was  ended. 
Two  different  time  periods  were  used  for  the  tests,  one  hour  and 
about  twenty  minutes.  The  reference  temperature  of  the  unused  room 
during  these  short  time  periods  was  again  swept  through  values  between 
the  outside  temperature  and  the  temperature  of  the  used  rooms . The 
results  of  these  runs  are  displayed  graphically  in  Appendix  F.  As 
the  graphs  show,  the  results  are  much  the  same  for  the  non-steady  state 
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case  as  for  the  steady  state  case.  The  curves  are  less  linear  bvit 
the  tread  is  upward  with  increasing  reference  temperature.  The  graph 
for  the  unused  time  of  twenty  three  minutes  shows  that  the  rooms 
do  not  have  time  to  cool  off  during  this  short  time  period,  so 
the  amount  of  heat  saved  is  very  little. 

The  proposal  of  finding  a median  temperature  which  would  be  used 
for  the  reference  temperature  for  the  unused  rooms  in  order  to  Affect 
a decreased  heat  usage  when  compared  with  other  temperature  settings 
for  these  rooms  has  been  thoroughly  investigated.  A graph  of  the 
heat  use  equation  developed  for  the  simplified  model  will  illustrate 
this  temperature . Mere  we  assume  heating  and  cooling  to  show  the 
opt imum  curves . 

TOTAL  HEAT  USED  VERSUS  UNUSED  ROOM  TEMPERATURE 


1 


■ — — T—l — —■ 


UNUSED  ROOM  TEMPERATURE 
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The  most  energy  was  saved  by  setting  the  unused  room  temperature 
to  some  value  below  the  steady  steady  state  temperature  in  that  room. 
This  chapter  has  given  some  analytical  reasons  why  the  controller 
was  programmed  to  work  the  way  it  was.  It  has  been  shown  that  the 
best  temperature  setting  for  a room  that  i6  not  to  be  used  for  a 
period  of  time  greater  than  an  hour  is  somewhere  below  the  steady 
state  temperature  ror  that  room.  If  the  contents  of  the  room  make 
it  necessary  to  maintain  that  room  above  the  steady  state  temperature 
then  the  lowest  setting  possible  is  the  most  conservative  of  energy. 
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CHA1  f LR  7 
CONCLUSION 

In  Chapter  6 certain  conclusions  were  reached  concerning  the 
temperature  setting  of  an  unused  room.  This  temperature  setting 
was  determined  to  be  any  temperature  lower  than  the  steady  state 
temperature.  It  other  constraints  enter  the  problem,  such  as  the 
need  to  keep  the  contents  of  the  room  above  the  freezing  level  or 
at  some  other  temperature  which  the  user  desires,  the  lowest 
temperature  consistent  with  these  constraints  is  the  most  energy 
efficient  for  the  unused  room.  Since  these  constraints  as  well  as 
the  decision  about  which  rooms  will  be  used  and  the  time  periods 
of  this  use  must  be  determined  by  the  person  or  persons  who  live 
in  the  house,  a model  cannot  be  made  of  the  use  of  rooms.  However, 
by  taking  what  are  typical  work,  sleep,  and  other  activity  patterns 
and  making  several  use  charts,  some  conclusions  may  be  reached 
concerning  the  amount  of  energy  saved  by  these  use  habits.  Two 
tvpical  use  charts  were  made(see  Appendix  G)  and  used  to  write 
programs  for  the  controller. 

The  first  use  chart  is  made  to  represent  the  use  habits  of  a 
working  couple  whose  children  attend  school  all  day.  The  house  will 
be  unused  from  the  time  the  family  all  leave  until  they  return  in 
the  afternoon.  The  bedrooms  are  not  used  until  the  evening  but 
the  other  rooms  are.  After  retiring,  the  living  room  and  kitchen 
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are  not  used  again  until  the  early  morning.  This  use  chart  can  be 
programmed  with  four  time  intervals  and  the  appropriate  temperature 
settings  for  each  room.  It  is  assumed  that  the  lowest  temperature 
to  which  an  unused  room  can  cool  is  35  degrees  F.  The  rooms  that 
are  being  used  are  set  to  75  degrees  F , and  the  outside  temperature 
is  at  20  degrees  F.  This  set  of  parameters  was  programmed  on  the 
H-8  and  used  to  control  the  house  simulation  for  an  equivalent  real 
time  of  24  hours.  The  total  heat  used  is  then  compared  with  the 
heat  used  to  heat  the  house  for  24  hours  with  all  rooms  being  heated 
to  75  degrees  F for  the  entire  period,  also  with  an  outside  temperature 
of  20  degrees  F.  For  this  use  chart  a savings  of  48.6  % was  realized. 
This  represents  a substantial  saving  of  both  energy  and  money  used 
to  purchase  that  energy.  The  model's  inherently  incorrect  results 
are  effectively  cancelled  by  examining  the  net  percentage  saving. 


The  second  set  of  use  habits  are  somewhat  less  saving  of  energy. 
The  case  of  having  someone  who  does  not  leave  the  house  is  considered 
in  this  use  chart.  At  night  the  bedrooms  and  the  bathroom  are  heated 
and  the  rest  are  unheated.  During  the  day  the  opposite  is  true, 
all  rooms  except  the  bedrooms  are  heated.  When  programmed  and 
compared  in  the  same  manner  as  the  first  use  chart,  a savings 
of  31.9  X was  realized.  These  tests  show  that  the  controller  does 
perform  its  task  in  a spectacular  way.  For  a house  such  as  the 
one  modeled  here  the  money  saved  by  a 30  X reduction  in  heating 
costs  spread  over  a winter  such  as  that  in  1976-1977  would  be 
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quite  a substantial  sum.  This  alone  would  suggest  that  the  digital 
controller  be  adopted  by  anyone  whose  house  is  configured  in  such 
a way(all  electric,  each  room  heated  separately)  that  the  installation 
costs  are  low.  The  builder  of  a new  home  would  also  be  able  to 
use  the  controller,  the  fixed  costs  of  installation  being  cut 
down  by  the  inclusion  of  the  controller  in  the  original  house 
during  construction.  There  are  other  ways  that  the  controller 
and  the  simulation  can  be  used. 

The  practice  of  lowering  the  setting  of  thermostats  was  one 
means  of  lowering  the  use  of  energy  in  homes.  How  efficient  was 
this  practice?  To  find  out,  the  simulation  was  set  up  in  a manner 
similar  to  the  tests  run  for  the  use  charts.  In  this  case  the 
rooms  were  all  heated  to  the  same  temperature  during  the  entire 
period.  When  the  amount  of  heat  used  over  a 24  hour  period  with 
the  thermostats  set  at  65  degrees  F is  compared  to  the  reference 
heat  use,  that  of  75  degrees  over  a 24  hour  period,  a savings  of 
17.3  7.  is  realized.  The  10  degree  difference  is  a little  extreme 
but  the  test  shows  that  a significant  amount  of  energy  is  saved. 
Another  application  of  the  simulation  is  the  measuring  of  the 
heat  saving  affected  by  adding  insulation.  To  change  the  insulation 
in  the  simulation,  the  governing  equations  are  changed  to  reflect 
the  new  thermal  transmission  coefficients  in  the  walls.  The 
changes  are  reflected  in  the  analog  simulation  in  the  changing 
of  potentiometer  settings.  If  a large  building  were  being  simulated, 
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a computer  could  be  programmed  to  change  the  potentiometer  settings 
with  an  input  of  insulation  values.  A program  which  takes  as  an 
input  the  values  of  the  insulation,  computes  new  equation  coefficients, 
computes  new  potentiometer  settings  and  then  sets  the  potentiometers 
could  be  easily  written  and  implemented  on  the  digital  portion  of 
a hybrid  computer.  A parameter  sweep  involving  the  input  of  a 
minimum  and  maximum  value  of  insulation  for  walls,  for  ceilings, 
for  roofs,  and  for  windows  could  be  written.  This  program  would 
check  the  heat  used  by  every  combination  of  insulation  values  and 
reach  what  would  be  an  optimum  insulation  configuration . Common 
sense  tells  us  that  this  would  be  the  thickest  insulation  everywhere, 
but  perhaps  this  is  not  the  case.  The  interior  walls  for  example 
could  be  either  totally  insulated  or  uninsulated  without  changing 
the  interface  to  the  outside  at  all.  The  study  of  the  most  efficient 
way  to  heat  our  buildings  is  one  of  the  most  important  ways  to  help 
stop  the  Impending  energy  crisis.  The  heating  of  living  spaces  is 
the  one  single  major  use  of  heat  that  cannot  be  curtailed.  People 
must  have  a livable  environment  or  they  will  freeze  to  death.  The 
use  of  digital  controllers  in  homes  will  save  some  of  the  energy 
now  used  to  heat  homes  and  the  use  of  the  house  simulation  could 
be  used  to  devise  insulation  and  construction  methods  to  save 
even  more  energy. 

Yet  another  advantage  to  using  the  digital  controller  is  the 
possibility  of  using  the  controller  in  a time-share  mode  to  control 
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other  things  in  the  house.  Since  the  rooms  need  to  be  sampled  only 
about  once  every  several  minutes  and  the  decision  to  turn  the  heater 
on  or  off  takes  approxiamtely  a few  microseconds  the  microprocessor 
is  idle  for  a large  percentage  of  the  time.  This  idle  time  could 
be  used  to  control  such  tilings  as  the  lights  in  the  house,  the 
burglar  and  fire  alarms,  the  operation  of  various  appliances, 
and  many  other  specialized  tasks.  The  use  of  the  controller  to 
control  the  peak  power  load  of  small  industry  is  another  possibility. 

A program  which  monitored  the  electrical  load  and  made  decisions 
about  whether  a machine  could  be  turned  on  at  that  time  or  at  a 
later  time  could  keep  the  peak  power  level  below  the  rate  cut-off 
and  save  the  company  money.  A central  microcomputer  for  a house 
which  controlled  the  temperature  in  individual  rooms,  the  running 
of  appliances,  and  many  other  tasks  in  the  home  could  be  the  trend 
of  the  future.  As  the  microprocessor  industry  becomes  more  well 
developed  and  the  prices  go  down,  this  universal  computer  control 
concept  ma«v  become  a realltv.  The  dreams  of  the  earlv  computer 
pioneers  have  been  surpassed  again  and  the  machines  which  were 

once  the  province  of  the  scientist  are  now  the  tools  of  the  common  j 

I 
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man.  Such  uses  as  have  been  investigated  here,  the  control  of  the 
temperature  in  a house  for  energy  conservation,  will  become  commonplace 


in  years  to  come. 
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This  Appendix  contains  the  uMes,  i 1 lustrat ions , and  lists  of 
equations  referred  to  in  Chapter  The  house  wttich  was  simulated 
is  drawn  with  side  and  end  views,  the  materials  used  in  the  construction 
of  this  house  are  listed  and  the  respective  lumped  heat  transfer 
coefficients  are  assigned  to  the  materials.  A complete  list  ot 
the  simulation  equations  is  also  included,  with  each  space  being 
examined  in  detail  with  respect  to  the  heat  flow  and  temperature 
relationships  to  one  another  and  to  the  outside  temperature. 


kA 


dT 

3X 


V/he  re 


tt  - heat  transfer  oer  unit  ti  -:«  ( t) ; 
ct 

A - awa  of  the  section  through  which 
heat  l 3 flowing; 

dT  - temperature  difference  causing  heat  fl 

dX  - length  of  path  through  which  heat  '.'low 
measured  in  the  direction  of  heat  "low 
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a proportional i ty  factor  called  therrrvi 
cor.ductivi  tv . 


room  A 


room  i 


room  k 


room  !i 


room  l"1 


room  H 


room 
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Exterior  walls 


u “ .24  BTU/hr- f t -deg  F 
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Lumped  Heat  Transfer  Coefficients 


W»hkI  Sujinn 
(CUpl»o*kr\l) 

■Ml  *W  W*K 
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. y i; ; J 
Ujpiisl* 


pane  112,  t.i 
po ! • it  ion  5 .c  . 


Int  erior  walls 


u - .33  BTU/hr- ft  -deg  F 
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I n,,,,  I'  1 


page  15h,  table  ■♦-?> 
parti  t ion  5 . . 


Floors 


u - . 58  BTU/hr- ft  -deg  F 


1§i  - ’ 


page  15S,  table  *-11 
floor  2 .C . 


Cei linv 


.12  BTU/hr- ft~-deg  F 


C'7-.  . _ • ' 

: C*  ■»“ 


page  157,  table  4-10 
ceiling  5.F.. 


u - .55  BTU/hr- ft2-deg  F 


l 

\ 


Inside  Doors  u - .64  BTU/hr- ft  -deg  F 
(hits ide  Doors  u - .48  BTU/hr- ft2 -deg  F 


Windows  u - 1.243  BTU/hr- ft  -deg  F 


page  161,  table  4-14 
roof  l . I . 


page  163 


page  163 
page  164 


lB . H.  Jennings,  Environmental  Engineering,  Analysts  and  Practice  (New 
York:  International  Textbook  Company, l d701 , pp . 152-164. 
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ROOM  A 


INTERFACE 

AREA( ft2) 

COEFFICIENT 

A to  Outside 

232  (1  door,  3 windows) 

78.294 

A to  E 

29.333  (1  door) 

15.57 

A to  F 

82.666  (1  door) 

33.17 

A to  H 

48.  (1  door) 

21.73 

A to  B 

72 

23.76 

A to  Floor 

176 

102.08 

A to  Ceiling 

176 

21.12 

dT/dt  (deg.  F/min)  = -.1909  * Ta  + .0506  * Tout  + .0101  * Te  + 
.0124  * Tf  + .014  * Th  + .0153  * Tb  + .0659  * Tfl  + .0136  * Tct 
+ .0006455  * Heater  A 


1 
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ROOM  B 


INTERFACE 

AREA (ft2) 

COEFFICIENT 

R to  Outside 

192  (1  door,  2 windows) 

62.676 

B to  A 

72 

23.76 

B to  H 

48  (1  door) 

21.73 

B to  I 

72 

23.76 

B to  Floor 

135 

78.3 

B to  Ceiling 

135 

16.2 

f i 

h 

! 

DT/dt  (deg.F/min)  « -.1905  * Tb  + .0527  * Tout  + .02  * (Ta  + Ti) 


r 


h 
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ROOM  C 


INTERFACE AREA(ft2) COEFFICIENT 


c 

to  Outside 

176  (2  windows) 

54.276 

c 

to  E 

53.33 

17.6 

I j c 

to  F 

58.66 

19.36 

4 

c 

to  H 

24  (1  door) 

13.81 

c 

. i 

c 

to  K 

88  (1  door) 

34.93 

to  Floor 

143 

82.94 

hi 

, C 

j 

to  Ceiling 

143 

17.16 

I 


! 

, 

j 

dT/dt  (deg.  F/min)  = -.1907  * Tc  + .014  * Te  + .0153  * Tf  + 
.011  * Th  + .0278  * Tk  + .0431  * Tout  + .0659  * Tfl  + .0136  * 
Tci  + .0007945  * Heater  C 
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ROOM  E 


INTERFACE 

AREA (ft2) 

COEFFICIENT 

E to  Outside 

24 

5.76 

E to  C 

53.333 

17.6 

E to  A 

29.333  (1  door) 

15.57 

E to  Floor 

11 

6.38 

E to  Celling 

11 

1.32 

dT/dt  (deg.  F/rnin)  = -.4816  * Te  + .1818  * Tc  + .1608  * Ta 


.0595  * Tout  + .0659  * Tfl  + .0136  * Tci 
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INTERPACK 
F to  A 
F to  C 
F to  H 
F to  FI oor 
F to  Ceiling 


ROOM  F 


AREA  1 1 1 ~ ) COEFFICIENT 


82.666  (l  doorl 

33.17 

‘'8 . non 

14.16 

24  (1  door 1 

13.81 

22 

12.76 

•j't 

2.64 

dT/dt  (deg.  F/mlnl 


.4221  * Tf  + 


.1713  * Ta  + .10  * Tc  + 


0713  * Th  + . 0659  * Tf 1 + .0136  * Tci 


5 


1 
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ROOM  H 


INTERFACE 

AREA(ft2) 

COEFFICIENT 

H to  A 

48 

(1 

door) 

21.73 

H to  F 

24 

(1 

door) 

13.81 

H to  C 

24 

(1 

door) 

13.81 

H to  K 

24 

(1 

door) 

13.81 

H to  D 

32 

(1 

door) 

16.45 

H to  I 

56 

(1 

door) 

24.73 

H to  B 

48 

(1 

door) 

21.73 

H to  Floor 

48 

27.84 

H to  Ceiling 

48 

5.76 

dT/dt  (deg.  F/min)  = -.3771  * Th  + .05143  * (Ta  + Tb)  + 

.0327  * (Tf  + Tc  + Tk)  + .0389  * Td  + .0577  * Ti  + .0659  * Tfl 
+ .0136  * Tci 
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ROOM  I 


INTERFACE 

AREA ( ft2) 

COEFFICIENT 

I to  Outside 

56  (1  window) 

19.458 

I to  B 

72 

23.76 

I to  H 

56  (1  door) 

24.37 

I to  D 

72 

23.76 

I to  Floor 

63 

36.54 

I to  Ceiling 

63 

7.56 

53 


T’  X13C» 


■rswmi*. 
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ROOM  K 


INTERFACE 

AREA (ft2) 

COEFFICIENT 

K to  Outside 

24 

5.76 

K to  D 

88  (1  door) 

34.93 

K to  H 

24  (1  door) 

13.81 

K to  C 

88  (1  door) 

34.93 

K to  Floor 

33 

19.14 

K to  Ceiling 

33 

3.96 

dT/dt  (deg.  F/min)  = -.3874  * Tk  + .0198  * Tout  + .1203  * Td 
+ .0475  * Th  + .1203  * Tc  + .0659  * Tfl  + .0136  * Tci 


+ .002583  * Tt  + .00562  * Td  + .001353  * Tk  4-  .001068  * Th 


+ .00902  * Tf  + .000451  * Te  + .005863  * Tc  + .1803  * Tout 


APPENDIX  A 


ATTIC 


AREA (ft  ) 


INTERFACE 


C0EEF1CT EN 1 


886.81  (roof I + 


Attic  to  Outside 


166.28  (walls'* 


527.652 


21.12 


Attic  to  R 


Attic  to  1 


Attic  to  D 


16.44 


Attic  to  K 


Attic  to  H 


17.16 


dT/dt  (dep.  F/mir0  - -.2118  * Tel  + .007215  * Ta  + .00553  * Th 


i 


ki 


[ 

\ ' 

[ I 
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APPENDIX  B 


This  Appendix  contains  a listing  of  the  general  purpose  digital 
simulation  program,  DIGISIM,  developed  by  E . E.  Mitchell.  This 
program  is  configured  to  simulate  the  model  of  a house  used  in 
this  project  and  is  set  up  to  print  out  a list  of  numbers  corresponding 
to  the  temperature  in  rooms  A and  C and  the  total  heat  used  over 
the  time  period,  not  specified  in  the  listing.  This  program  was 
of  great  use  in  removing  the  errors  in  the  model,  since  the  digital 
simulation  integration  technique  is  sensitive  to  small  errors  in 
the  model  and  these  small  errors  caused  the  simulation  to  "blow  up." 

By  running  the  model  again  and  again  until  the  temperatures  were 
maintained  at  a steady  value  with  no  heater  inputs,  the  errors  were 
removed  . 


I 


i 


i 


I 


IBOli  COPY  FURNISHED  TO  DD.C  ^ 


i'  i g i s :.  > 


D1GISIM  MG  CAL  SIMULATOR  — F.K.  MI  "VHF1.1 


a**  REMEMBER  — RUN  FORTRAN 

s 

* LINKS  400-650:  INITIAL  TIFF  DA"' A FLOCK 

* L.  I ' 3 750-1000:  SYSTEM  DYNAMIC  EQUATIONS 

* -ON  ••  S 1000-1500:  OU”  PUT  STATEMENTS 

- L : 1700-1850:  TERMINAL.  CALCULATION  REGION 

x 

' B A : • I v DATA:  TMAX  , DEL,"' , NSAMPL 

" Pa.  . STATEMENTS:  X ENTG  RL  ( XD , X ( 0 ) ) , CALL  OUTPUT  X , "X"  ) 

* A!  ' NAL  ’ ’ ONS  : STEP  , FRSTOR  , SON  DOR' , PULSE  , SQUARE  , PK! . A Y , D IGTA1  , ' ' 'T  " 

P"’RA  I N 

" FO:  ;•  HE  DETAILS  CONTINUE  LISTING 

* 

*#**  INITIAL  TIME  REGION 

* 

* This  PROGRAM  HAS  AN  INITIAL  TJMK  DATA  REGION  EXTENDING  FROM  LINES  400-7' l 

* TUTS  SPACE  IS  USED  TO  DEFINE  CONSTANTS,  INITIAL  VALUES  AND  OTHER  NON 

* CHANGING  PARAMETERS  REQUIRED  BY  YOUR  PROGRAM.  ALSO  INITIAL  CALCULATIONS 

* MAY  PF  MADE  HERE  (IF  YOU  DEFINE  THE  REQUIRED  rARAME”ERS). 

* IN  ADDITION,  SPECIAL,  C0NS'nANTS  FOR  BUILT  IN  DIG  I ST  M FUNCTIONS  ARE  DP  IN!' 

* HERE.  FOR  INSTANCE , COMMONLY  DEFINED  TERMS  T.N  '"HIS  AREA  ARE: 

* HE AD=  ” ” ONE  LINE  OF  80  OR  LESS  CHARACTERS  BETWEEN  THE  QIUV'E'T 

* THIS  IS  TYPED  OUT  AS  A HEADING  FOR  the  RUN 

* DEL?  INTEGRATION  STEP  SIZE  ( DEI. "TO . 05  DEFAULT' ) 

* ""i-'iU  TOTAL  SIMULATED  RUN  TIME  ( "'MA X - I DFKAUL"' ' 

* N.  '."T,L  OUTPUT  SPACING-  ANSWERS  PRINTED  EVERY  D EL"’  * NS  A OP  I 

* TIME  INTERVAL  (NSAMPL-?  DEFAULT) 

* ZK  D.U,  "’AL  CONTROLLER  GAIN 

* 7.A(  1-4)  ,ZB(  1-4)  MORTAL  CONTROLLER  COEFFICIENTS  ( Z A ( ) . 7,  B ( ) - 0 DK  • ' 


"■'SAMP 

TLAG 

TO 

TKPI.O" 


PPL ANE 


DIGITAL  CONTROLLER  SAMP  TIME  (TSAMP  =0  DLLFAin.’" ) 

DELAY  FI  D FOR  TIME  I’ El. AY  (TLAG-0  DEFAULT) 

INI"'1  AL  TIMT-'  TK  CHER  whaN  0 (TO-O  DFFA'H."' ) 

IF  TR  PLOT--  1 '"HE  OUTPUT  RESPONSE  IS  WHl'CEN  TN'"0  A F !!.'•' 
SAVED  BY  YOU.  IT  IS  IN  '"11F  PROPER  FORM  FOR  n.OTTING 
ON  THE  TKKTRONIC  GRAPHIC  TERMINALS .( '"KP1.0T  :0  DFEAUir) 
>■  P PLANE-:  1 THE  "’KPLO"'  FILE  IS  WRITTEN  TO  MAKE  A PDAS'- 
PLAN'  Pi.O"'  0F  THE  FIRS WO  CAL'.  OU'"PUt  VART  API. KS 


"»*  DYNAMIC  PEG  ION 

* 

* THL  PROG  HA  I HAS  A DYNAMIC  REGION  FROM  STATEMENT  NUMBERS  750-1000.  IN  ”'H  I s 

* AREA  thf,  SYS’" KM  DYNAMIC  , CONSISTING  OF  DIFFERENTIAL  AND  ALCFBRA"TC 

* EQUA'"rONS,  A R F DEFINED. 

* ' N HE  DYNAMIC  REGION,  THE  FUI.L  REGULAR  FORTRAN  LIBRARY  IS  AVAIL. 1 ABLE. 

* FOR  "HAT  MA'"  I'.  1'"  IS  AVAI.T  ABLE  ANYWHERE  IN  '"up  PROGRAM.  IN  ADD  I"'  1 ON 

* "HE  •OLLOWtKC  ONS  A R h FURNISHED:  ( X - I N B L: " , V r OU  P t'  ’" ) 


TIME 


r 


ii 


lVn^cTV 

*»«&>** 
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DTGISIM  (continued) 


= ENTGRL  (X , XO)  — XO=  INITIAL  VALUE  - DECIMAL  POIN"  REQUIRED  IF  NUMBER 
= STEP  ( T 1 ) — Y=  0 FOR  T<T1,  Y=1  FOR  ">=T1 
= PULSE ("  1 ,T2)  — Y = 0 EXCEP"  Y=  1 WHEN  " 1 < = "< r"V 

= SQUARE  ( T 1 , T2 , T3  , P ) --  SQUARE  WAVE  WITH  AMPLITUDE  OF  Y=1.  LEADING  EDGE 

OF  FIRS"'  PULSE  A’"  PULSE  WID"F  IS  "'1  "0  "2, 

REPEATS  EVERY  P SECONDS. 

= FRSTOR ( X , A , B , C , D ) — FIRST  ORDER  TRANSFER  FUNCTION  DEFINED  AS 

(A»S+B)/(C*S+D)  C IS  NO"’  ZERO 

= SCNDOR ( X,A,B,C,D,E,F)  — SECOND  ORDER  TRANSFER  FUNCTION  DEFINED  AS 

(A*S*S+B*S+C)/(D*S*S+E*S+F)  D IS  NOT  ZERO 
= DELAY  ( X , TL AG ) — TIME  DELAY,  Y = X AFTER  X IS  DELAYED  "'LAG  SECONDS 
= DIGTAL ( X ) — DIGTAL  FILTER  OR  CONTROLLER  SIMULATOR . DEL"  SHOULD  PF 

<=TSAM/10.  TSAM--SAMPLE  INPUT-OU"PU"  TIMES  OF  "HE  DIG"AL 
CONTROLLER.  TRANSFER  FUNCTION  IS: 

(ZK»(ZA(  1)*X+ZA(2)»X-1  + ZA(3)*X-3)  - ZB(2)*Y-1 
-ZB ( 3 ) * Y-2  )/ZB(1) 

X= CURRENT  INPUT, X-1  = LAST  INPU",ETC.  DI"TO  FOR  Y 
= LIMIT ( U 1 , X , U2 ) — Y=U1  IF  X<U1,  Y=U2  IF  X>U2,  Y = X OTHERWISE 

PTRAIN(TI)  — Y=  1 IF  T = N*71;  THIS  IS  A PULSE  "RAIN,  UNIT  PULSE 
EVERY  N * T 1 SECONDS 


***  OUTPUT  REGION  : 

THE  PROGRAM  HAS  AN  OUTPUT  REGION  FROM  S" AXEMEN"  NUMBER  1200-1500.  IN  "HIS 
REGION  THE  VARIABLES  "'0  BE  PRINTED  AND  PLOTTED  ARE  DEFINED.  TO  OUTPUT 
X,  XDOT  AND  X2DOT , THE  STATEMENT'S  ARE: 

i 

CALL  OUTPU"(  X,  ,rX"  ) 

CALL  OUTPUT ( XDOT , "XDO"" ) 

CALL  OUTPUT ( X2D0T , "ACCEL " ) 

ANY  8 CHARACTERS  MAY  BE  USED  BETWEEN  THE  " ”,  THESE  A h THE  COLUMN 
HEADINGS  OVER  "HE  PRINTED  VALUES. 

* NOTE  **  TIME  IS  OUTPUT  AUTOMATICALLY 

* NOTE  **  THE  FIRS"  6 CALL  OUTPU"  VARIABLES  ARE  PLO""ED 

***  TERMINAL  CALCULATION  REGION 

THE  PROGRAM  HAS  A TERMINAL  REGION  FROM  STATEMENT  NUMBERS  1700-1850.  THIS 
AREA  IS  USED  TO  MODIFY  PARAMETERS , EC".  A"  THE  END  OF  A RUN.  ONE  CAN 
THEN  GO  BACK  AND  .1 N "EG RATE  AGAIN  AND  AGAIN.  A "YPIGAL  PROBLEM  "HA" 

WOULD  USE  THIS  AREA  IS  A POUND  ARY  VALUE  PROPI.FM. 


IF  YOUR  PROGRAM  CONTAINS  THE  S"ATEMEN'r 
GO  TO  9000 

CONTROL  IS  TRANSFEREE  BACK  TO  THE  TOP  OF  "HE  INITIAL  AREA 
IN  A LIKE  MANNER,  THE  ST ATKMFNT 
GO  TO  9100 


MiniriwuiyfrT^Mi 


-sssss^2. 


mCnc^ 


D 1 G I S I !■  i ( c o 1 1 ' m i u c d ) 

* TAKES  CON '."HOI.  'r0  THE  BOTTOM  OK  "'HE  INTTTAI  AREA,  ("'OP  OK  DYNAMIC  AREA) 

# 

***»  END  OF  INSTRUCTIONS 

* * » 


I M P L I C I'1'  REAL  ( I-N  ) 

I \ "EGER  ICZ  , NICZ  , I TZ,NI"7. , TTMAX  ,KFRR  , T AZ  , NZZ  , KPZ  , NSAMKL  , K"*Z  , T 7 Z , JJZ  ,T7 
1 EC ! ' E 1 C N ,r  , N V A R , IDIGL. , ICN  , I ADCNT  , T.SZZ 

COMMON  7. 0 ll  7 , X V Z ( 6 , 4 0 ? ) , I C N , K V A R , T , ,n  0 , T M A X , D L ' L 1 , U D E 1 , 1 , I C Z , N 1 C Z , I T Z , N T - ' 
& t'rMAX, KERR, TAZ(  100'  ,Y1PZ(  100)  ,Y4PZ(  100)  , X 1 P Z ( 1 00 ) , X4  FZ  ( 100)  ,N  77  ,EEZ  , 

& NEAMDL ,KTZ 

COMMON  /liOLD/IADCN"  , ITQA  , JCN  " , IQCNT , HLD  ( BOO  ) 

C 0 MMON  /OU"’ P2  / I I 7, , J J 7. , A Z Z , I SZ Z , A M X Z ( ?r> ) , AMN  Z ( 25 ) , H F AD , N P LOT  ,*rKP L0’~  , 
&PKLANE 

COMMON/DIG  ITA/.IDIGL  , ZA  ( 4 ) , ZB  ( 4 ) , TSAM  , ZK  , NCN,p 
COMMON /'TEMP  1 /SWA  , SUB  , SWC  , SWD  , SWH  , SKI 
CO M MO N / T F M P2/TA  , "’B  , TC  , TI) , TH  , "'I  , TOU"' 

CHARACTER  PDQ,PDQ2 , IPQ , AZZ ( 25) , HEAD»80 

REAL  LIMIT 

INTEGER  TKPLOT 

DATA  ICN.IIZ, JJZ/0,0,0/ 

ZOUTrICNTrlDICL- 0 
”NPLOT=0 
N PLOT =6 
DO  100  17.  = 1.  A 
Z.  A ( T Z ) = 0 . 

ioo  zn ( i z ) =0 . 

ZK-O. 

T ADCN’prO 
DO  108  17- 1 , EOO 
108  HLD(IZ)-0 
1 DO  ? X Z , 1 100 
? I A Z ( I.  'x - 0 
OOOGCOXTINL'C 
T 0 - 0 . 

DELT= .05 
NICZs  10 
N'  "7. -If 

***'  INITIAL  REGION  ***  LINES  400-650 


■'ll  ' 

jl 1 

1 


jHlSPAGfE 


fiW.lSBBUW1100  
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DIGIS.LM  (continued) 


* 

* 

«»*« 

IF ( TKPLOT  .EQ.  0)GG  TO  9100 
WRITE( ZOUT ,67 ) 

87  FORMATC  WHAT  IS  YOUR  SAVED  FILE  NAME") 

READ  ( ZOUT  , 1 4 ) IPQ 
OPENFILE  3, IPQ 
REWIND  3 
ENDFILE  3 
9 100KTZ=0 
T = TO 

hDELT=0. 5*DELT 

KERR=0 

ICZ-  1 

IThAXsNICZ 

4 ITZrl 

5 NZZ-0 

****  DYNAMIC  REGION  #**  LIMES  750-1000 

TEDr-,48l6*TE+.  1818*TC+.  1608*TA+. 0595*T0UT+ . 0659*TFL+ . 0 1 36*TCI 


TFD=-. 422 1 #TF+ . 1713*TA+. 1 *TC+ . 07 1 3*TH+ . 0659*TFL+ . 0 1 36*TCI 

TKD  = -.  3874*TK+.0196*T0UT+.  1 203*TD+ . 0475*TH+ . 1203*7C+.0659*TFL+.01 36*  .’Cl 


TAD=-.  1909*TA+.05U6*T0UT+.0101*TE+.0214*TF+.014*7H+.015  3*TP+.0659*TFL 
&+ . 0 1 36*TCI+QINA 

TBDr-.  1905*TB+.0527*T0UT+.02»TA+.02*TI+.0183*TH+  .06!j9*T?L+.0136*TCI 
A+QINB 

TCD  = - . 1 907*TC+ . 0 1 4 *TE+ . 0 1 5 3*TF+ . 0 1 1 * TH+ . 027C  * TK+ . 04  3 1 *T0UT+  . 0659*TFL 
A + .0136*TCI  + C?INC 


TDD-- . 1 900*TD+ . 029*  TK+ . 0 1 97* TI+ . 0 1 36*TH+ . 0482* TOUT+ . 0659 *TFL 
&+.0136*TCI+QIND 


TID  = -.2443*TI+.04  29*TB+.04  4*TH+.0428*TD+.0351  M’OUT+.0659*TFL+.01 36*TCI 
&+QINI 

7CID=-. 2 1 16*TCI+. 00721 5*TA+.00553*TB+. 00258 3*Ti+. 00562* TD 
&+. 001 353*7K+. 00 1 968* TH+.000902*TF+. 00045 1*TE+. 005863*TC+. 180 3*TOUT 

THD=-. 3771 *TH+ .0514 3*TA  + . 05 1 4 3*TB+ . 0 327*TF+ . 0327* TC+ . 0 J27*TK+ . 0 389* TP 
&4.0577*TI+.0659*TFL+.0136*TCI+QINH 


TCI=ENTGRL ( TCID , TC.IOO) 


IUIS  PACE  IS  BEST  QUALITY  PRACTICABLE 
FROM  COPY  FURNISHED  TO  DDC  


- DIG! SIM  (continued) 

TH=ENTGRL  ( '."HD  , ”'HOO ) 

TP=ENTGRL ( TBD , TBOO  ) 
TlrEN7C.RL('7ID,7I00) 

TF  = ENTGRL ( TFD , TFOO) 

TA  - EN  TG  RL  ( D , T A GO ) 

TI  rEN",GKL(7ED,TEOO) 
TP=ENTGRL(TDD ,7D00) 
Tk-ENTGRL("KD  ,'°KOO) 

TC'  - EN^GFL  ( mCD  , 7CU0 ) 

CALL  THERMS 

i QIN’S=SWB*HEAT*.8'»16 

i.' i N i = S W I * H E A 7 * 1.80  35 
Q.INA  = SWA*HFA7*,64  55 
QIND=SWD*HEAT».8293 
0INC  = SWC*HEAT».79i*5 

Q'"A  = ENTGRL  ( QINA  ,0.0) 
QTB=EN?GRL ( QINB , 0.0) 
QTI=ENTGRL(QINI ,0.0) 
QTC=ENTGRL(0INC,0.0) 
CTD-ENTGRL ( QIND , 0 . 0 ) 
GTDD-QINn+QINA+QINC+QIND+QINT 
Q^-'ENTCRL  (QTDD  ,0.0) 


Kb AT = 18 

TOUT=TFL=TFMP2 
TEMP? =20 

TC I00=TB00=7I 00  = TFG  O = TA00  = TE00  = TD00  = TK00  = 
ATHOOsTEMPI 
TCOOsfjO 
TEMPI =75 

IF (KERR  .NF.  0)  GO  TO  9 
ITZ  = ITZ+ 1 

IF  ( ITZ  . I.E  . ITMAX ) GO  TO  5 
KPZsO 

***  OUTPUT  REGION  **  LINES  1200-1500 
CALL  OUTPUT ( QT  , "TOT AL  HFA'1’"  ) 

CALL  OUTPU~(TC , "TC"  ) 

CALL  OUTPU,r,(TA  , "TA"  ) 

JJZ=JJZ+2 

IF(ICZ  .EQ.  0 ) GO  TO  6 
ICZ  = 0 

ITMAX = NITZ 

6 IF (T  .L~.  (7MAX+.00001))G0  TO  68 
IF ( TKPL07  . NE . 0)G0  TO  7 
WRI7E(Z0U~,  1?) 

12  FORMAT ( " WAN"’  TO  INCREASE  "’MAX"  ) 


jj»S«SS5£S& 
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PIGIStM  (continued) 


READ(ZOUT, 14)IPQ 
14  FORMAT ( V) 

I F ( I P Q .NE.  "YES")  GO  TO  7 
WRI TE( ZOUT , 13) 

13  FORMAT ("  TMAX=") 

READ (ZOUT, 14 )TMAX 
GO  TO  6 
68  KTZ  = KTZ+ 1 
T=KTZ 

T=T*DELT+TO 
GO  TO  4 

9 WRITE ( ZOU'r  ,10) 

10  FORMAT ( 5X  , " INTEGRATION  ASSIGNMENT  ERROR") 
STOP 

7 ISZZ  = 2 

***  TERMINAL  REGION  ***  LINES  1700-1850 


IF(TKPLOT  .EC.  0)CALL  RANGE 
I F ( TKPLOT  ,EQ.  0)  CALL  PL0T2 
IF(TKPLOT  .NE.  0)  CALL  TEKPL" 

STOP 

END 

FUNCTION  ENTGRL ( X , YO ) 

COMMON  0UT,XVZ(6 ,402)  , ICN’t  ,NVAR  , T , "0  , t^X  , DELm , HDFLt  , IC,NIC,  IT, NX”, 
&ITMAX,  KERR,  IA(  100)  ,Y1(  100)  ,Y2(  100)  , X 1 ( 100)  tX2(  1 00 ) , N , KP  , NSAMPL  , KT 
N =N+  1 

IF(IC.EQ.O)  GO  TO  20 
IF(IT.NE.I)  GO  TO  10 
IF(IA(N).NE.O)  GO  TO  50 
IA ( N) = 1 
GO  TO  15 

10  IF( IA (N ) . EQ . 0 ) GO  TO  50 
15  Y2N= YO 
GO  TO  40 

20  IF(IA(N).EQ.O)  GO  TO  50 
IF(IT.NE.I)  GO  TO  30 
1 PREDICTOR 

Y1N=Y2(N) 

X IN- X2( N) 

Y1(N)=Y1N 
X1(N) rXIN 
Y2N=Y1N+X1N*DELT 
‘ CORRECTOR 

30  Y2N=Y1(N)+(X1(N)+X) *HDELT 
40  Y2 ( N ) = Y2N 


ms  PAQE  IS  BEST  QUALITY  PRACTICABLE 
mm  cott  FwtajdtD  toddq  -- 


PIGI3IM  (continued) 

X2 ( N) =X 
EN'rGRL  = Y2N 
RETURN 
r>0  KF.RR=  1 
RE’-UHN 
END 

FUNCTION  S"‘EP(T1) 

COMMON  OUT,XVZ(6,402) , ICN™ ,NVAR ,T 
Yr  1. 

IFC~.LT.T1 ) Y-0. 

STEP* Y 
UPTURN 
END 

FUNCTION  PULSE (T1,T2) 

COMMON  OUT , XVZ ( G , 402 ) , ICNT , NVAR , T 
Y-0. 

IF(T.GE.T1.AND.T.LE.T2)  Y= 1 . 

PUI.SEsY 

RETURN 

END 

FUNCTION  SQUARE(Tt ,T?,T3) 

COMMON  OUT  , X VZ  ( 6 , 4 Oc? ) , I CN”' , NVAR  , T , n'0  , '"MA X , nFL”’ , HPFL"’ 

Y=0, 

TK= AMOD ( T-T 1 , T3 ) 

IF (TM  .GE.  0.  .AND.  TM  .LE.  T2)  Y = 1 . 

SQUARErY 

RETURN 

FND 

FUNCTION  F R STO R ( X , A , B , C , D ) 

COMMON  /FRST/  °7.  ,Z 
l F ( C . EQ . ZERO ) GO  TO  1 
SZ=X-D*Z/C 
Z = ENTGRL (SZ , 0 . ) 

FRSTORr ( A *SZ  + B*Z )/C 
RETURN 

1 PRINT  2 

2 FORMAT  ( 5X  , "FIRS'”  ORDFR  TRANSFER  FUNCTION  Oir  OF  ORDER”) 
STOP 

END 

FUNCTION  SC N DOR (X,A,B,C,D,E,F) 

COMMON  /SCND/  S2Z,ST.  ,Z 
I F ( D . F.Q . ZERO ) GO  TO  1 
S2ZnX-(H#SZ+F*Z)/D 
SZ  = ENTGRL (S2Z ,0.  ) 

Zr ENTGRL ( SZ , 0 . ) 

SCNDOR= ( A*S2Z+B*SZ+C*Z )/D 
RETURN 

1 PRINT  2 

2 FORMAT ( GX , "SECOND  ORDER  ’’’RANSFER  FUNCTION  0Um  OF  ORDFR”  ) 
STOP 


! 
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D LG  S ■ *'  (continued) 

END 

SU  'ROUTINI-  OUTPUT  ( X , A 1 ) 

CH*.:  ACTER  * 1 ,A(25) ,HFAD*80 

COD.  MO N OIT  , XVZ  ( 6 , <40 7 ) , ICNT  , NVAR  ,T  , '"0  , TMAX  , DEI  , HDEL'r  , TC  , NI C , TT  , N I - , 
ATT  /AX  , KER:  , 7 A ( 100  ) , Y 1 ( 100)  , Y,?(  100  ) , XI  ( 100  ) , X2{  Ur  ) , N , KP  , NSAMP!  , KT 
COMMON  /0UTP.?/IC,JJIA,IS:'Z,AMXZ(25)IAMNZ(?S),HJ-A  .NPLOT.TKFLO  ,!  LANF 
CO!'  MCN  /OU'-pl/  V(25)  , IP  , K.'.SE  ,KS,TSAV  , ICN 

* 

I F ( T . N E . '"O’-  -0  TO  10 
KS  = 0 

I F ( J J .GV  0)  GO  <10 

11  = 11+  t 

A v I T. ) = A 1 

l S o L = 1 

NVAR  = Ml N ( II , 6 , NFLOT) 

GO  TO  <40 

10  IF(KP  .NE.  0)  GO  TO  <10 
I F ( K S .EC.  0)  GO  TO  60 
00  KS=KS+1 

I F ( K S .GE.  NS AM PL  ) KS  = 0 
<10  KP=  KP+  1 

IF ( KS  .NE.  0)  RETURN 
TSAV  = T 
V ( KP)  = X 
RETURN 


R I N 

T LAST  T 

IME 

PERIOD 

CON 

TINUE 

[ r v 

'T-  = ICNT+  1 

DO 

1 20  1=1, 

KVA 

R 

XV 

(I, ICNT) 

= 

vci) 

I F ( 

TKPL.OT  . 

NE. 

0 }(■■■'  ''"0 

IF (ICN  .EC. 

0) 

WR I ' E(OU 

IF  ( 

ICN  .EQ. 

0) 

prin" 

if; 

ICN  .1*". 

91 

) GO  V 7 

ICN  = 0 

WR.ITE(OUT,  9 15) 

WR.  "EC OUT  ,Q00) 

70  CONTINUE 

kt  ■ ~ r:i  n c o,n ) 

I F ( J J .PC.  2 ) WRITE ( OUT  ,910) (A(I) ,1= 1 ,KT1 ) 
7 F ( 1 1 .GT.  5) WRITE (OUT, 905) 
V;RT’-E(0UT,90'.')lCfr,-'AV  , (V(T)  ,1=1  ,K"1  ) 

ict;=icfi+ 1 

IF  (II  .LE.  5)  GO  ,r0  100 

► 

ICN  = ICN  + 1 
SO  KTisKTl  + 1 

KT,  - MIN(m  + ;f  ,11) 

IF(JJ  . EQ.  ?)WPITE(0UT,9^0)  ( A ( I ) ,I  = in 

WR  - TE(OU-,9J»0)(V(I)  , T = K "r  1 .K’V) 


> K ) 
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DIGIS1M  (continued) 

ICN=ICN+ 1 
KT1  = KT2 

IF (II  .GT.  K T2)  GO  TO  80 
TOO  JJ=4 

CALL  RANCH 
GO  TO  <20 

90''  FORMAT ( 1 H 1 ,//////) 

90  > FORMA T( 1H  ) 

910  FORMAT ( / 1 1 X , "TIME"  , 3X  , 5 ( ?X  , A8  , 2X ) ) 
918  FORMAT (////) 

920  FORMAT ( 1 H , I 3 , 1 X , 6G 1 2 . 5 ) 

9 30  FORMAT ( / f 1 8X , 5 ( 2X , Al , 2X ) ) 
ylJO  FORMAT ( 1 H , 15X.5G12.5) 

END 

SUBROUTINE  PL0T2 


N=  NUMBER  OF  POINTS 
NVARr  NUMBER  OF  PLOTS 


V=  NVAR  X N MATRIX;  EACH  COLUMN  CORRESPONDS  TO  A VECTOR 
COMPOSED  OF  POINTS  TO  BE  PLOTTED 


MAX  NVAR  = 6 

MAX  N = ttOO  WITHOUT  CHANGING  DIMENSIONS 


COMMON  OUT , V ( 6 , U02 ) , N , N VAR 

COMMON  /OUri,.Wll  , JJ , A2 ,ISZZ , AMXZ(29)  , AMNZ(28)  , HEAD  , NPLOT  . TKPLO"' , 
&PPLANE 

DIMENSION  SF ( 7 ) , IS ( 7 ) , H ( 7 ) , IC ( 7 ) 

CHARACTER  BLANK , DOT , S7R , MINS , HEAD*80 
CHARACTER  X(8),LINE(  120) ,A2(>5) 

BLANK="  " 

DOT=" I " 

STR="+" 

M J NS= 

X(  1 )="♦■• 

X ( 2 ) = " * " 

DO  9999  1=1, 7, 2 
9999X ( I ) =X( 1 ) 

DO  9998  1 = 2, 8, 2 
9998X(I)=X(2) 


* 

» 


! 
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DIG1SIM  (continued) 


DO  1 1=1 , N V A K 
H( I ) =0 
SF ( I ) =0 . 0 
IS(I)=I 
i ic (i) si 

DO  6 I = 1 , N V A R 
V ( 1 , 1 ) = V ( 1 , 1 ) + 1 . E-6 
DO  6 J = 1 , N 
IF(V(I,J)  ) 3 , 3 , ^ 

3 ICQ)  =2 

I K ( A BS  ( V ( I , J ) ) ~H  ( 1 ) ) 6 , 6 , S 

5 H ( I ) = ABS( V ( I , J ) ) 

6 CONTINUE 


GENERATE  SCALE  FACTORS 


DO  7 1= 1 , NVAR 
A = 60/ ( NV A R* IC ( 1 ) ) 

SF(I)sA/(H(I)  + . 1 * H ( I ) ) 

A = ( 30  * ( IC ( I ) - 1 ) +60  * ( I- 1 ) ) /NVAR 
DO  7 Jsl.N 

7 V(  l, J)=5F(I)»V(If J)+A 


WRITE  READING 


W R I TE ( OUT , SO 0 ) N V A A 
IF ( NVAR-** ) 8 , 8 , 9 
8 NS=  1 
NQ=  1 
N P = N V A R 
GO  TO  10 
0 NS  = 2 
NQ=  1 
NT  = 14 

10  DO  11  K=  1 , NS 

WRITE( OUT , S0 1 ) ( A2 ( I ) , I = NQ , NP) 
WR1TF,(  OUT  , S02 ) ( SF(  ! ) , 1-NQ  , NT  ) 
NQ  = S 

11  MPsNVAR 
WRITE(OUT,50S) 

I’Rl  NT , HFAD 

WR T TE  ( OU"’ , SO  3 ) 


SB 


i 


V 


J 


tf 

ft 

ft 


MARK  AXIS 
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IISIM  (continued) 


LLNTH  = 01 
K = 0 

DO  12  1=1  .LLNTH 
LINE(I)=MINS 

WHITE(OUT,504)K,  (L1NE(I)  ,1=  1, LLNTH) 


BLANK  THE  LINE 


DO  13  1=1  , LLNTH 
LINE(  I ) = BLANK 


PLOT  THE  VARIABLES 


DO  15  J=  1 ,N 

LINE ( 1 )=DOT 

DO  14  1=1 , N V A R 

JK= ( 60  *I)/NVAR  +1 

JP= JK-( 30  *(IC(I)-1)  )/NVAR 

JL=V(I,J)+2.5 

LINE( JK ) =DOT 

LINE ( JP ) =00'” 

LINE( JL) =X( 1 ) 

WhITE(0U7, 504)  J,  (L.INE(K)  ,K=  1 , LLN,rH  ) 
DO  15  K= 1 , LLNTH 
LINE( K) =BLANK 


RETURN 


500  FORMAT  ( ////5X  , "SYSTEM  PLO,r  IN "/5X  , IP  , 5X  , " V A Rl  ABL.FS"/ ) 

501  FORMAT ( BX  , " V ARIABL.ES;  " , 5X , 4 ( A8 , 6X ) ) 

•Or  FORMA"* ( 3X  , '•SCALE  FACTORS"  , 4 (G  1 4 . 5 ) ) 
j03  F0RMAT(///1X,6( 10X," 1")) 

504  FORMAT ( 1 X , 1 4 , 6X  , 1 2 1 A 1 ) 

•05  FORMAT (// ) 

END 

FUNCTION  DELAY ( C , TL AG ) 

CHARACTER  A ( 25  ) , HE Al) *80 

COMMON  OUT  , XVZ  (6,402)  , 1CNT  ,NVAR  ,T  ,T0  ,TMAX  , DEL"* , HDEL"* , IC  , NIC  , 

A I T t NIT ,ITMAX , KERR , l A ( 100) , Y 1 ( 100),Y4(  100 ) , X 1 ( 100),X4(  100),N,KP, 
ANSAMPL ,KT 
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DIGISIM  (continued ) 

COMMON  /0UTP2/  U , J J , A , ISZZ  , AMX  ( 25 ) , AMN  ( 25  ) , HEAD  , NPLO'" , '”KPI  v"  . 
ftp PLANK 

COMMON  /0U,nP1/  V(?5)  ,IP,KASE  ,K  f?SAV,ICN 
COMMON  /HOLD/  I ADCN'r  , 1TQA  , JCNT  , IQCNT  ,HI.D(  500  ) 

IF(1ADCNT.GL.  1 ) GO  TO  102 

IADCNTb  1 

ITQA-  1 

JCNT=  0 

IQCNT=0 

10?  IF(T.GT.TLAG)  GO  TO  103 
DLL AY  = HLD ( ITQA ) 

HLD(ITQA) =C 
IQCNT = IQCHT+ 1 
GO  "'0  10  A 

103  IP (JCNT. GE. 1 ) GO  TO  100 
ITQA=  1 
JCNT=  1 
GO  TO  107 

106  IF( ITQA. LT. IQCNT)  GO  TO  107 
J C N T = 0 

107  DELAY=HLD( ITQA ) 

HLD ( ITQA )=C 

10 '1  ITQA  = ITQA+  1 . 

RETURN 

END 

FUNCTION  D.T.GTAL  ( U ) 

COMMON  OUT ,XVZ  (6  , iJ02)  , ICN'1' , N VAR  , T , TO  , TMA  X , DEL.’1' , HDELT  , I C , NIC  , I T , NT~ 
ftlTMAX  , KERR  , I A ( 1 00  ) , Y 1 PZ  ( 101))  , Y4PZ(  1 00 ) , X 1 PZ  ( 100 ) , XU  PZ  ( 100 1 , N , K P , , K 
ANSAMPL 

COMMON/DIG  IT  A /IDIGL  , 7,  A ( 4 ) ,ZB(4)  ,TSAM  ,ZK  , NCNT 
SAVE  Y , Y 1 ,Y ?.  Ill  , U n , U 3 » Y 3 


IF (IDIGL  .C".  0)  GO  TO  100 
IDIGL=  1 
NCNT-0 
U1  0. 

U2  = 0. 

Y 1 =0 . 

Y?--0. 

U3=0. 

Y3  = 0. 

* 

* ADJUST  DELT  TO  PE  INTEGRAL  MUL^IPLF  OF  TSAM 

* 

I F(  DELT  .cr.  TSAM/ 10.)  GO  'n0  50 
PP=1. 

20  RR  = TSAM/ PP/ DEL"' 

JJrRR 

I F ( ( RR-FLOAT  ( J J ) ) .L’r.  .001)  GO  TO  70 
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DILI  SIN  (continued) 

PPsPP+1 . 

1F(PP  . LE . 3.)  GO  TO  20 
PP=1.  ’ 

30  XDEL=?SAM/20./PP 

IF ( XDEL  .LT.  DELT)  GO  TO  60 
PP=PP+  1 . 

GO  TO  30 
60  XDEL=TSAM/ 10 . 

60  N S A M P L = N S A M PL*DE’LT/XD E L 
DELT=XDEL 
GO  TO  TOO 
70  DELTsDELT/PP 

NSAMPL=NSAMPL*PP 
100  CONTINUE 
* 

* HERE  FOR  NORMAL  OPERATION 

* 

I F ( APS  ( K'"*DELT  - NCN"'*TSAM ) .GT.  1.E-6  ) GO  TO  300 

IF ( IT  .EQ.  ITMAX ) NCN,n  = NCN,n+ 1 

Y-ZK*(ZA( 1 )*U+ZA(2)*U1+ZA(3)*U2+ZA(4)*U3) 

Ys ( Y-ZB ( 2 ) * Y 1 -ZB ( 3 ) * Y2-ZB ( 4 ) *Y3 ) /ZB ( 1 ) 

IF (IT  .LT.  ITMAX)  GO  TO  300 

U?  = U1 

U1  = U 

Y2  = Y 1 

Y1=Y 

* 

300  DIGTALsY 
RETURN 
END 

SUBROUTINE  RANGE 
CHARACTER  A (25) ,HEAD*80 
COMMON  OUT,XVZ(6 ,402) , ICNT 

Common  /OUTP2/  II  , J J , A , IS  , AMX ( 25 ) , AMN  ( 25  ) , HEAD  , NPLO'”  , TK'  i.OT, 
& PPL AN E 

COMMON  /OUTP 1 / V(?5) ,IP , K ASE , KS , TSA V ,ICN 

* 

IF(ICN,n  .GT.  1)  GO  TO  20 
DO  10  1=1,25 
AMX ( I ) = -1.E+6 
10  AMN ( I ) = 1.E6 
20  GO  TO  (30 ,60), IS 
30  DO  50  1=1,11 

AMX(I)  = AMAX1 (AMX (I ) ,V(I)  ) 

50  AMN(I)  = AMIN  1 ( AMN (I),V(I)) 

70  RETURN 

* 

80  WRITE (OUT, 900) 

DO  100  .1=  1 , II 

100  WHITE ( OUT , 9 10 ) I , A (T ) , AMX ( I ) , AMN ( I ) 


-1 
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DIGISIM  (continued) 


GO  TO  70 

900  FORMAT ( 1 HO  , " MAX  AND  MIN  VALUES") 

910  FORMAT ( 1 U ,I2,3X,A8,"  MAX  r ",G12.5,"  MIN  = ",012. 5) 
END 

REAL  FUNCTION  LIMIT  ( I1 1 , X , U2  ) 

LIMIT  = X 

IF(X  .GT.  U2 ) LI MI Tr U2 
I F ( X .LT.  U1)  LIMIT  = U 1 
RETURN 

END 

SUBROUTINE  ""  E K P L T 

* j 

* TH‘S  SUBROUTINE  PREPARS  A FILE  FOR  PLOTTING 

* ON  THE  TEKTRONICS  GRAPHIC  TERMINALS. 

* YOU  MUST  OPEN  A FILE  IN  YOUR  CATALOG 


CHARACTER  A ( 25 ) , HEA D*80 

COMMON  OUT,  X(  6, 402)  , ICNT  , NVAR  , T , TO  , TMAX  , DEL'1' , HDEL"’ , IC  , NTC  , IT  , N T ' 
&.ITMAX  , KERR  , I A ( 1 00  ) , Y 1 PZ  ( 100 ) , Y4  PZ  ( 100  ) , X 1 PZ  ( 100  ) , X4  PZ  ( 100  ) , N . KP  , 
&NSAMPL 

COMMON  /OU,rP2/  II,JJ,A,IS,AX(25)#AZ(25),HEAD,NPL0?,TKFL0T,r»M  v 


IF( PPLANE  .NE.  0)G0  ”0  70 
DO  50  I = 1 , N V A R 
K = 0 
T = TQ 

CO  20  J = 1 , ICNT 
WRITE(  3 , 500  ) T , X ( I , .J ) 

K = K + 1 
T r K*NSAMPL 
T = T*DELT  + TO 
20  CONTINUE 

WRITE( 3 » 5 10 ) 

50  CONTINUE 
RETURN 

70  DO  80  Jsl.TCNT 

80  WRITE( 2 » 500) X(  1 , J ) , X ( 2 , J ) 

WRITE( 3,510) 

RETURN 

500  FORMAT (El 2. 4 , El 2. 4) 

510  FORMAT ( " 1.E37  , 1 . E 3 7 ") 

END 

REAL  FUNCTION  PTRAIN(TI) 

COMMON  Z0Um,XVZ(6 ,402) ,ICNT,NVAR,T,TO,TMAX 
SAVE  P,J 
DATA  J/0/ 

PTRAINzO. 

P=J 


jRjuoonFvsausH®™01' 


mClia3l* 
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P = P*T 1 

IF (T  .GT . P)J=J+1 

I F(ABS(T-P)  .LT.  .0001 )P"RAIN= 1 . 

RETURN 

END 

SUBROUTINE  THERM 1 

COEMON/TEMP 1 /SWA , SHE  , SWC  , SWD  , SWH  , SWl 
C0EM0N/TEKP2/TA  .:TB  , TC  , TD  , TH  , TI  , TOUT 
R E F = 7 5 
LOW=35 

As. 5* ( REF-TOUT )+TOUT 
IK ( A. GT. LOW)  GO  TO  70 
TFLO=LOW 
GO  TO  71 

IF (A . GT.REF ) GO  TO  75 
TFLOsA 
GO  TO  71 
?FLO= REF 

IF(T.GT.240)  GO  TO  1 
GO  TO  3 

IF  VT.GT.720)  GO  TO  2 
RC  = RD  = RH=TFLO 
RA  = RI  = REF 
KBs LOW 
GO  TO  4 

IF(T.GT. 1200)  GO  TO  3 

RHsTFLO 

KC  = RD  = RI  = REF 

RA=RB=LOW 

GO  TO  4 

RA=RB=RH=TFLO 

RIeREF 

KC  = RD  = LOW 

CONTINUE 

IF ( TH .GT . RH ) GO  TO  845 
SWH  = 1 
GO  TO  846 
SWH  = 0 

IF(TB.GT.RB)  GO  TO  849 
SWB=  1 
GO  TO  850 
SWB  = 0 

CF( TI .GT .HI)  GO  TO  853 

SWIrl 

GO  TO  854 

SWIsO 

IF(TA.GT.RA)  GO  TO  857 
SWA  s 1 
GO  TO  858 
SWA  = 0 


* 


I 
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DIGISIM  (continued) 

858  IF(TD.GT.RD)  GO  TO  861 
SWD=  1 
GO  TO  862 

861  SWD  = 0 

862  1F(TC.G"\FC)  GO  TO  866 
SWC  = 1 

GO  TO  866 
865  SWC -0 
86b  CONTINUE 
RETURN 
END 

SUBROUTINE  THERM? 

COMMON/TtT.F  1 /SWA  , SWB  , SWC  , SWD  , SWH  , SWT. 
COMMON/ TEMP? /TA  , TB , TC  , TD  , TH  , '"l  , "’OUT 

REF=75 

RA=RB=RD=RH=Rl=REF 
i RC=50 

846  IF(TB.GT.RB)  GO  TO  849 
SWB=1 
GO  TO  850 
649  SWB=0 

850  IF(TI.GT.RI)  GO  TO  853 
SWI  = 1 
GO  TO  854 

853  SWI=0 

854  IF(TA.GT.RA)  GO  TO  857 
SWArl 

GO  TO  858 

857  SWAsO 

858  IF(TD.GT.RD)  GO  TO  86? 

SWD  = 1 

GO  "'0  86? 

861  SWD  = 0 

862  IF(TC.GT.RC)  GO  TO  865 
SWC  = 1 

GO  TO  866 

865  SWC=0 

866  CONTINUE 
RETURN 
END 


75 


F 


h 
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APPENDIX  C 

This  Appendix  contains  the  analog  computer  flow  charts  used  in 
programming  the  mathematical  model  of  the  house  on  the  EAI-681 
analog  computer.  The  flow  charts  have  been  presented  in  an  equation 
by  equation  basis,  because  of  the  great  number  of  components  involved. 
The  input  signals  into  the  summers  on  each  page  are  the  outputs  of  the 
integrators  on  other  pages.  The  disjointed  nature  of  these  diagrams 
make  interpretation  difficult,  but  the  maze  of  connecting  lines  that 
would  result  from  a connected  display  would  be  even  more  undecipherable. 
Also  included  in  this  Appendix  is  the  digital  portion  of  the  hybrid 
simulation.  The  various  programs  which  were  used  to  set  up,  control 
and  record  the  results  of  the  analog  simulation.  The  programs  included 
are  HEAT,  with  its  associated  subprograms,  PARA,  THERM,  TRAN,  and  0PT1. 
The  potentiometer  setting  program  P0T1  is  included  also.  These 
programs  were  all  developed  as  a direct  result  of  this  project. 

HEAT  is  the  main  program  which  calls  up  the  various  subprograms  in 
order  to  run  the  simulation.  PARA  sets  up  the  parameters  for  a run 
of  the  simulation,  THERM  is  the  heater  control  subprogram,  TRAN  is 
the  subprogram  which  transfers  data  to  and  from  the  analog  computer, 
and  OPTl  is  the  digital  controller  subprogram. 
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SUBROUTINE  PAJ? A < K ) 

INTEGER  OPTC 20 ) 

L)I  MENS  ION  JT  ( 20 ) ,P(b),KT(20  ) , 10(20),  IT(20) 
C()MM()N/mLKT/KT/BLKO/IQ/BLKR/JT/BI.K()/()PT 
COMMON/BLKP/P/BLK I /IT 
CALI.  HINITCIE) 


GO  TO  (bOO, 1000)  , K 
CONI'  INUE 


I 

I ' 

4 


WRITE (4, 100) 

REAO( 4 , 1 10)  LPT 

100  FORMAT ( 22H  NORMAL  PARAMETERS?  YES=1;  N0=0) 

110  FORM  Aid  2) 

IFILPT.EQ. 1 ) GO  TO  1 70 

WR I TEC  4.  120) 

READ( 4 , I TO ) J T ( 12) 

120  FORMAT ( 2bH  INPUT  THERMOSTAT  SETTING) 

130  FORM AT ( 1 5 ) 


WRITE (4, 140) 

READ (4 , 1 30)  PC  2) 

140  FORMATC  25H  INPUT  HEATER  SIZE  RTU/HR ) 


; 


i 


WRITEC4, 160) 

READ (4,1 30 ) KT  ( 2 ) 

160  FORMAT ( 3bH  INPUT  RUN  TIME  IN  SIMULATION-HOURS) 
KT( 1 1 )=JOO+KT( 2 ) 

WR IFF ( 4 , 16b) 

READ ( 4 , 1 30)  KT ( 10) 

16b  FORM AT(25H  INPUT  INITIAL  CONDITIONS) 

WR  ITE ( 4 , 167) 

READ ( 4 , 1 6B ) IOC  lb ) 

REAOC  4,168)  10(16) 

167  FORMAT ( 2 OH  INPUT  I Q( I b ) , IOC  1 6 ) ) 

168  FORMATC  2 I b ) 

GO  TO  300 


L l— 
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PARA  (cont.) 


170  KT( 1 1 )=7500 
KT( ! 0 )=7200 
P ( 2) =200 

JT(  1 2)=JT(8)=JT(4)=JT( 1 0 ) =JT( 2 ) =JT ( 3 ) =7500 
IQ (5 )= IQ ( 6 )=2000 
KT( 1 2) =0020 


300  CALL  SAM0( 7,  IE ) 

CALL  SAMOC 7, IE) 

CALL  SPOT (002 ,0050, 0002,  IF*) 
CALL  S POT ( 012,0100, 0002, IE) 


CALL  SPOT(010,KT( 11 ), 0002, IE) 

- 

i 

1 

U GO  TO  44 


1000  CONTINUE 


CALL  SP()T(  006,  JT(  20),  0002,1  E) 


WR  ITE (4,44) 

44  F()RMAT(  1 7H  ANALOG  OPERATING) 


RETURN 

END 
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SUBROUTINE  THERM 

JIM::  NS  ION  M( 5 ) , 1 5(  10),  IT(20  ) , I Q(  20  ) , P(  5 > , JT(  20 ) 
COMMON/BLKI /IT/BLKO/IO/BLKR/JT/BLKP/P 
R=  .5 
V.  ( I ) =4 

:/( 2)  = io 

r.U  3)  =3 
M(4) =8 
.4(5)  =2 
DO  5 L=l ,5 

IF ( IT ( M( L ) ) . GT. JT( M( L ) ) ) GO  TO  1 

IS  (M  ( L ) ) = 1 

GO  I'O  5 

IS (M( L ) ) =0 

CONTINUE 


1 0(4  )= I S ( 4 ) *P( 2 ) * 1 8. 035*R 
10(3 ) = 1 3 ( 1 0 ) *P( 2 ) *B. 4 1 6 

I0( I ) = (IS( 4 )*R  + IS( 10) +I5( 2)  +IS (8)  + IS( 3 ) )*P( 2) 
IQ ( 8 ) = IS ( 2 ) *P ( 2 )*7 . 945 
10(7  )=  13  ( 8 ) *P( 2 ) *8. 293 
10(2 )=IS(3)*P(2)*6. 455 
RETURN 
END 


SUBROUTINE  TRAN ( L ) 

DIMENSION  IT(20> ,10(20 >tKT(  20) 
COMMON /i<LK  I /IT/BLKO/IO/BLKr/KT 
CALL  H IN  IK  ID 


GO  iO  ( I ,3  ),L 
I DO  2 K=l , 10 


2 CALL  CRACS(K-1 ,IT(K>,IE> 

CALL  CRACS(0I2,KT( I ) , IE) 
CALL  CRACS (011,10(10), IE) 
GO  l'O  S 

3 DO  4 K=l  ,8 

CALL  LTDAS(K-1 , IQ( K) , IE) 

4 CALL  TLDA 

5 RETURN 

END 
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SUB  ROUT  I Nil  OPT  I ( K ) 

INTEGER  OPT (20) 

1)1, ML-.  NS  ION  IQ(  20  ) , KT(  20 ) , JT(  20) , IO(  20) 
COMMON  /BLK  T/KT/BLKR/JT/BLKO/  IQ/BLKO/OPT 
HI =7 S 00 
L()h=3500 

GO  10  ( 7 , b ) , K 
CONI l NUE 

I P ( K 1( I ) .or. 1200)  GO  TO  8 
JT(4  ) = JT  ( 2 ) = J T ( 8 ) =L()W 
JT( I 0 ) =JT ( 4 )=H I 


IF(KT(  I ) .Gr. 2700)  GO  TO  1 
GO  TO  3 

I F ( ;<  r ( I ) .G1'.4S00)  GO  TO  2 
JT(S )=LOU 

JT (3 )-jr ( 2 ) =JT( 4 ) =JT ( I 0 )=HI 
C^C ) i O 4 

CONI  I NUT 

J T ( 2 )= JT  ( 8 ) = JT  ( 4 ) =H  I 
JT  (3  )= JT  ( 1 0 )=L()W 
GO  i'O  4 

J I' ( 3 ) =JI(  I 0 ) = J T ( 4 ) =JT ( 2 )=JT(  8)  =1.0 w 
J T ( 3 ) = J I ( I 0 ) = JT ( 4 ) =H  I 
CONTINUE 
GO  TO  8 


J r ( U = J i’(l  0 ) = JT  ( 3 ) = JT ( 8 ) = jr ( 2 ) =7800 


RETURN 


I 
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POTl 


11=0002 


CALL  HI  NIK  IE) 

•CALL  SAM0(7,IE) 

CALL  SPOT (057, 005 1 , IT, I E) 
K=057 

IF(IE.NE.O)  WRITEM, I )IE,K 
CALL  SPOT! 023,0033, IT,  IE) 
K=023 

IF(IE.HE.O)  WR ITE( 4, 1 ) IE, K 
CALL  SPOT! 0)6, 0033, IT, IE) 
K=01  6 

IF(IE.NE.O)  WRITE! 4, I ) IE,K 
CALL  SP0T(020,0043, IT, IE) 
K=020 

IFIIE.NE.O)  WRITEM,  ))  IE, K 
CALL  SPOT(063,0!04, IT, IE) 
K=063 

IF(IE.NE.O)  WRITEM, I )IE,K 
CALL  SPOT (06 1 ,0022, IT, IE) 
K=06  I 

IF ( I E.NE.O  ) WRITE! 4, I ) IE,K 
CALL  SPOT (064,0) 43, IT, IE) 
K=064 

IF(IE.NE.O)  WRITE(4, 1 ) IE,K 
CALL  SPOT (062,001 5, IT, IE) 

IK=062 

IF(IE.NE.O)  WR ITE( 4 ,1 ) IE, K 
CALL  SP0T(0J2,0136, IT, IE) 
K=032 

IF(IE.NE.O)  WRITE! 4, I ) IE,K 
CALL  SPOT (060,1 900, IT, IE) 
K=060 

IF(IE.NE.O)  WRITE ( 4, 1 ) IE,K 
CALL  SP0f(072,0066, IT, IE) 
K=07  2 

IF (I E.NE.O)  WR ITE (4,1 ) IE, K 
CALL  SPOK  069,0050,  IT,  IE) 
K=069 

IF (I E.NE.O)  WRITE! 4, I ) IE, K 


I 


i 

I 

| 


H 


h 


■ 

. 
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POTl  (cont.) 

C SET  POTS  FOR  ROOM  B 

CALL  SP()r(  1 I 1 ,0020,  IT,  IE  ) 

K=  1 1 1 

IF(IE.NE.O)  WRITEC4, I )IE,K 
CALL  SPOr<095,0020, IT, IE) 
K=OvS 

I F ( IE.NE.O  ) WRITE! 4, 1 ) IE,K 
CALL  SPOT (OR I ,0183, IT, IE) 
:<=ori 

I E ( IE.NE.O)  WRITE(4, I )IE,K 
CALL  SPOI'C  110,1  V0<3,  IT,  IE) 

K=  1 1 0 

IFC IE.NE.O)  WRITEC 4, 1 )IE,K 
CALL  SP()T<  I 15,0527,  IT,  IE) 

K=  I I 5 

IFC IE.NE.O)  WRITEC 4, 1 ) IE,K 
C SET  POTS  FOR  ROOM  E 

CALL  SPOT (018,0181 , IT, IE) 
K=0\  8 

IF (IE.NE.O ) WRITE (4, 1 ) IE,K 
CALL  SP()i'(  101,0160,  IT,  IE) 
K=1 0 1 

IFC IE.NE.O)  WRITE ( 4, ) )IE,K 
CALL  SPOT( 100,4821 , IT, IE) 

K=  1 00 

IFCIE.NE.O)  WRITEC4, 1 ) IE,K 
CALL  SP()T(  1 02,0595,  IT,  IE) 
K=102 

IF(IE.NE.O)  WRITE(4, 1 ) IE,K 
C SET  POTS  FOR  ROOM  K 

CALL  SP0TC017, 1 203, IT, IE) 
K=01  7 

IF ( IE.NE.O)  WRITEC4,  I)  IE.fC 
CALL  SP()T(033,0120,  IT,  IE) 
K=033 

IF( IE.NE.O  ) WR ITE( 4, 1 ) IE. K 
CALL  SPOT (037 ,0047, IT, IE) 
K=03  7 

IF ( I E.NE.O)  WR  ITE ( 4 , 1 ) 15, K 
CALL  SP0T(030,3870, IT, IE) 
K=030 

IF ( IE.NE.O  ) WRITEC  4, 1 ) IE,K 
CALL  SPOT (031 ,01 98, IT, IE) 
K=03  1 

IF( IE.NE.O)  WRITEC4, 1 ) IE,K 
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POTl  (cont.) 

c SHI'  POTS  FOR  ROOM  F 

CALL  SPOT  (073, 01  71  , IT, IE) 
K~07  3 

I F( I E . NH . 0 ) WR ITE ( 4 , I ) IE,K 
CALL  SP0r(013,0l00,IT,IE) 
K=01  3 

IF( I E.NH.O)  WRITEC4, I )IE,K 
CALL  SPOT ( 070 ,42 1 9, IT, IE) 
K=070 

IF(IE.NE.O)  WR I TEC  4, I ) IE,K 
CALL  SP0T(090,0713, IT, IE) 
K=090 

IF( I E.N5.0)  WRITEC4, 1 )IE,K 
C SET  POTS  FOR  ROOM  Cl 

CALL  SPOT (055 ,0010, IT, IE) 
K=0S5 

IF(IE.NH.O)  WRITE (4, 1 ) IE,K 
CALL  SPOT( 047,0070, IT, IE) 
K=04  7 

IF(IE.NH.O)  WR ITE( 4,1  ) IE,K 
CALL  SP()T(  I 12,0045,  IT,  IE) 

K=  1 1 2 

I F ( I E . ME  . 0 ) UR ITE( 4, I) IE, K 
CALL  SPOT (031 ,001 3, IT, IE) 
K=08  I 

IF(TE.NE.O)  WRITE(4, 1 )IE,K 
CALL  SP0r( 083,0005, IT, IE) 
K=083 

IF(IE.NE.O)  WRITEC4, 1 )IE,K 
CALL  SPOT (027 ,0060, IT, IE) 
K=02  7 

IF( I E. ME .0 ) WR I TE ( 4, ! ) IE,K 
CALL  SPOT (054, 0060, IT, IE) 
K=054 

IF(IE.NE.O)  WR ITE ( 4, 1 ) IE, K 
CALL  SP()T(045,00I9,  IT,  IE) 
K=04  5 

IF(IE.NE.O)  WRITE(4, I )IE,K 
CALL  SPOT( 097,0026, IT, IE) 
K=0V7 

IF(IE.NE.O)  WRITE(4,  I ) IE,K 
CALL  SP0T(080,021I , IT, IE) 
K=0H0 

IF ( IE.NE.O)  WRITE(4, I ) IE,K 
CALL  SPOT ( I IP, 0180, IT, IE) 
K=H8 

IF( IE.NE.O)  WR ITE ( 4 , I )IE,K 
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POT1  (cont.) 

C SET  POi'S  H)R  ROOM  I) 

CALL  SP()r(  040,01  36,  IT,  IE) 
K=04  0 

IF(IE.NE.O)  WR ITE( 4, I ) IE,K 
CALL  spore  042,001 0,  IT, IE  ) 
K=04  2 

IF(  IE.NE.O  ) WRITE(  4,  1 ) IE,.< 
CALL  SP0T(087,0029, IT, IE) 
K=0H7 

IF( IE.NE.O)  NR  ITE(4 , 1 ) IE,K 
CALL  SP()f(085, 1 894,  IT,  IE) 
K=0b5 

I E ( I E . NE . 0 ) WR ITE (4,1  ) I E , K 
CALL  SPOT ( 079,0482, IT, IE) 
K=079 

IF (IE.NE.O)  WR  ITEM,  I ) IE,K 
L.  SET  POTS  FOR  ROOM  I 

CALL  SP0r(034,0351 , IT, IE) 
K=034 

IF  ( I E.N5.0)  NR  ITEM,  1 ) IE,K 
CALL  SP()T(035,0440,  IT,  IE) 
K=035 

I F ( IE.NE.O ) NR ITE ( 4 , 1 ) IE,K 
CALL  SP0T(067,0043, IT, IE) 
K=067 

I F( IE.NE.O)  WR ITE ( 4 , 1 ) IE, K 
CALL  SP()T(065,2447,  IT,  IE) 
K=065 

IF (IE.NE.O)  WRITEM,  I )TE,K 
C SET  POTS  FOR  ROOM  H 

CALL  SPOT(O46,003J, IT, IE) 
K=04  6 

IF(IE.NE.O)  WR ITE( 4, 1 )IE,K 
CALL  SP()T(07.I  ,0577,  IT,  IE) 
K=07  I 

IF(IE.NE.O)  WR ITE( 4, 1 ) IE,K 
CALL  SP0r(077,0389, IT, IE) 
K=077 

IF< IE.NE.O)  WR ITE (4,1  ) IE, K 
CALL  SPOT( 105,0052, IT, IE) 

K=  105 

I F(  IE.NE.O)  WR  ITEM,  1 )IE,K 
CALL  SP0T(075,3772,IT, IE) 
K=075 

IF(  IE.NE.O)  WRITEU,  1 ) IE,K 
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POTl  (cont.) 

C Str  POTS  FOR  ROOM  c 

CALI.  SP0r(052,0l 42,  IT,  IF  ) 
k=0S2 

IF(IE.NL'.O)  WRITH(4,  I )JE,K 
CALL  SPOl (022,01 12, TT, IF) 
K-022 

IF< IE.NE.O)  WRITE! 4, 1 ) IE, K 
CALL  SPOT (053,0028, IT, IF) 
K=OS  3 

IH< I IE.NE.O)  WR ITE ( 4 , I ) IF,K 
CALL  SP()1'(  07<s,0432,  IT,  IF) 

K ® 07  ft 

IF(IE.NF.O)  WRITFC4, I ) IF,K 
CALL  SPOT (OS  1, 00 15, IT, IF) 
K=OS  I 

1 F ( IE.NE.O)  WR ITF ( 4 , I ) IE,K 
CALL  SPOT (050,0192, IT, IE) 
K»0‘)0 

IF< IE.NE.O)  WRITE! 4, I ) IE,K 
I FORMAT (2X, I3,2X,I3) 

END 
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This  Appendix  contains  the  specification  sheets  for  the  interface 
unit  components.  These  specification  sheets  are  simplified  versions 
of  those  put  out  by  the  manufacturers  of  the  components.  They  are 
not  meant  for  use  by  someone  who  is  building  such  a device  but  for 
the  reader  to  gain  insight  into  the  function  and  use  of  the  interface 
device.  The  connection  of  the  components  as  well  as  the  Interaction 
between  them  is  more  easily  understood  when  a basic  knowledge  of 
the  components  is  gained.  The  schematic  for  the  interface  unit 
is  included  in  this  Appendix  also.  The  text  of  Chapter  4 contains 
a discussion  of  the  use  of  this  device. 


APPENDIX  D 
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SPECIFICATION  SHEET  1 

ADC-EK8B  Monolithic  Analog-to-Digital  Converter 
The  ADC-EK8B  Is  a low  power,  Integrating  ana  log- 1 o-d igi t a 1 
converter  fabricated  on  a single  monolithic  chip  using  CMOS  technology. 
The  circuit  employs  a charge  balancing  integrator,  current  switch, 
comparator,  clock  counter,  data  counter,  and  control  logic  circuitry 
to  implement  conversion.  The  charge  balancing  integration  technique 
gives  high  linearity  and  noise  immunity  along  with  inherent  monotonicit 
resulting  in  no  missing  codes.  Output  data  appears  in  parallel 
form  on  latched  outputs  which  are  CMOS,  low  power  TTL,  or  low 
power  Schottkv  TTL  compatible. 


PIN 

FUNCTION 

PIN 

FUNCTION 

1 

NC 

13 

Reference 

2 

NC 

14 

Analog  input 

3 

NC 

15 

Amp  1 1 f i e r out 

4 

NC 

16 

Zero  adjust 

5 

Bit 

0 

out  (MSB'' 

17 

Bias 

6 

Bit 

1 

out 

18 

-5  volt  power 

(in) 

7 

Bit 

2 

out 

19 

+5  volt  power 

(in) 

8 

Bit 

3 

out 

20 

Ground 

9 

Bit 

4 

out 

21 

Start  convert 

(in) 

10 

Bit 

5 

out 

22 

E . 0 C . (status ) 

11 

Bit 

6 

out 

23 

Data  valid 

12 

Bit 

7 

out 

24 

NC 

Rin  - FSR/10  A - lOv/10  A ■=  1M 
Rref  M vref/2°  A = +5V/20  A » 250  K 


2ADC-KK8B  Series 
Datel  Systems,  Inc. 
1020  Trunpike  Street 
Canton,  Mass  02021 


- 
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SPECIFICATION  SHEET  2 

CD4051A  COS/MOS  Analog  Multiplexer 

The  RCA  COS/MOS  analog  multiplexer,  CD405IA,  is  a digitally 

controlled  analog  switch  having  a low  "on"  impedance  and  a very 

low  "off"  leakage  current.  Control  of  analog  signals  up  to  fifteen 

volts  peak-to-peak  can  be  achieved  bv  digital  signal  amplitudes 

of  three  to  fifteen  volts.  The  multiplexer  circuits  dissipate 

extremely  low  quiescent  power  over  the  full  Vpp-Vpj.  supply  voltage 

range,  independent  of  the  logic  state  of  the  control  signals.  When 

a logic  "1"  is  present  at  the  "inhibit"  input  terminal  all  channels 

are  "off".  CD4051A  is  a single  8-channel  multiplexer  having  three 

binary  input  signals,  A,  B,  and  C,  and  an  "inhibit"  input.  The  three 

binary  input  signals  select  one  of  eight  channels  to  be  turned  "on" 

3 

and  connect  the  input  to  the  common  output . 


PIN  FUNCTION PIN  FUNCTION 


1 

Channel  4 (in’) 

9 

Control 

C 

(in) 

2 

Channel  6 (in') 

10 

Control 

B 

(in) 

3 

Common  output 

11 

Control 

A 

(in) 

4 

Channel  7 (in’) 

12 

Charnel 

3 

(in) 

5 

Channel  5 (in) 

13 

Channel 

0 

(in) 

6 

Inhibit 

14 

Channel 

1 

(in) 

7 

VEE 

15 

Channel 

2 

(in) 

8 

VSS 

16 

VDD 

Inhibit 

C 

B 

A 

"ON"  Channel 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

1 

0 

2 

0 

0 

1 

1 

3 

0 

1 

0 

0 

4 

0 

1 

0 

1 

5 

0 

1 

1 

0 

6 

0 

l 

1 

1 

7 

1 

- 

- 

- 

none 

3 RCA  Solid 

State 

1 Z4 

DATABOOK  Series 

COS/MOS 

Digital 

Integrated 

Circuits 

RCA  Solid  State, 

Box 

3200, 

Somerville,  N.  J.  0887b 
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SPECIFICATION  SHEET  3 
CD4010AE  COS/MOS  Hex  Buffer/Converter 
This  buf fer/convert er  provides  conversion  ranges  from  COS/MOS 
logic  operating  at  positive  three  volts  to  positive  fifteen  volts 
supply  level  to  DTL  or  TTL  logic  operating  at  positive  three  volts 
to  positive  six  volts  supply  level.  Conversion  to  logic  levels 
greater  than  positive  six  volts  is  permitted  providing  Vcc (DTI. /TTL) 

A 

is  greater  than  Vpp(COS/MOS) . 


PIN 

FUNCTION 

PIN 

FUNCTION 

1 

^CC 

9 

D (in) 

2 

A (out) 

10 

D (out) 

3 

A' (in) 

11 

E (in) 

4 

B (out) 

12 

E (out) 

5 

B (In) 

13 

NC 

6 

C (out) 

14 

F (in) 

7 

C (in) 

15 

F (out) 

8 

Ground 

16 

VDD 

4 


RCA  Solid  State  '74.  DATABOOK  Series 
COS/MOS  Digital  Integrated  Circuits 
RCA  Solid  State,  Box  3200,  Somerville,  N.  J.  08876 
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SPECIFICATION  SHEET  4 

High  Performance  Operational  Amplifier  pA74l 
The  ^A741  is  a high  performance  operational  amplifier  with  high 
open  loop  gain,  internal  compensation,  high  common  mode  range  and 
exceptional  temperature  stability.  The  J/A741  is  short-circuit 
protected  and  allows  for  nulling  of  offset  voltages. 


PIN  FUNCTION  PIN  FUNCTION  

Offset  null 
Output 
V+ 

NC 

^ TTL  Integrated  Circuits  Data  Hook 
Compiled  bv  0 . H.  Perry  III  LT,  USN 
Weapons  and  Systems  Engineer ing  Depavment , 

United  States  Naval  Academy 


Oftset  null 
Inverting  input 
Non- inverting  input 
V- 
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SPECIFICATION  SHEET  5 
Hex  Buffer/Driver  N7407 

The  54/7407  and  54/7417  Hex  buf fer/driver  features  standard 
TTL  inputs  with  non-inverted  high  voltage,  high  current  open 
collector  outputs  for  interface  with  MOS , lamps,  or  relavs.  The 
54/7407  maximum  output  is  thirty  volts  and  the  54/7417  maximum 
output  Is  fifteen  volts. ^ 


PIN 

FUNCTION 

PIN 

FUNCTION 

1 

A (in) 

8 

D(out) 

1 2 

1 

A (out) 

9 

D (in) 

3 

B (in) 

10 

E (out) 

-j  4 

B (out) 

11 

E (in) 

a 5 

C (in) 

12 

F (out) 

6 

C(out) 

13 

F (in) 

i 7 

I 

Ground 

14 

VCC 

) 

,(  ^TTL  Integrated  Circuits  Data  Book 

Compiled  by  0.  H.  Perry  III,  LT,  DSN 
Weapons  and  Systems  Engineering  Department 
Dnited  States  Naval  Academy 

I 

\ 
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MULTIPLEXER  SIGNAL  CONNECTION  SCHEMATIC 


r*\«jt  H PCtTM  l 
A*'1  O Si 
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EK8B-A1V  StCNAl.  CONNECTION  SCHEMA!'  !( 


Vvj.  K>-  Vs,,  .As.  ; \v*0'V  * r 

C «. t t «.  $ ft 


O ♦'  y 
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CONTROL  SIGNAL  CONNECTION  SCHEMATIC 


AW 
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PIN  DESIGNATIONS  FOR  INTERFACE  BOARD 


1. 

Control  bit  6 (out) 

A. 

-18  volts  supply 

2. 

Control  bit  5 (out) 

B. 

Analog  signal 

7 (in) 

3. 

Control  bit  4 (out) 

C. 

Analog  signal 

6 (in) 

4. 

Control  bit  3 (out) 

D. 

Analog  signal 

5 (in) 

5. 

Control  bit  2 (out) 

E. 

Analog  signal 

4 (in) 

6. 

Control  bit  1 (out) 

F. 

Analog  signal 

3 (in) 

7. 

Control  bit  0 (out) 

H. 

Analog  signal 

2 (in) 

8. 

+10  volts 

; supply 

J. 

Analog  signal 

1 (in) 

9. 

Control  bit  7 (out) 

K. 

Analog  signal 

0 (in) 

10. 

NC 

L. 

Control  bit  6 

(in) 

11. 

NC 

M. 

Control  bit  5 

(in) 

12. 

Control 

bit  7 (in) 

N. 

Control  bit  4 

(in) 

13. 

Control 

bit  0 (in) 

P. 

Control  bit  3 

(in) 

14. 

ADC  bit 

7 (out)  LSB 

R. 

Control  bit  2 

(in) 

15. 

ADC  bit 

6 (out) 

S. 

Control  bit  1 

(in) 

16. 

ADC  bit 

5 (out) 

T. 

ADC  data  valid  (out) 

17. 

ADC  bit 

4 (out) 

U. 

Multiplexer  control  A (in) 

18. 

ADC  bit 

3 (out) 

V. 

Multiplexer  control  B (in) 

19. 

ADC  bit 

2 (out) 

W. 

Multiplexer  control  C (in) 

20. 

ADC  bit 

1 (out) 

X. 

ADC  start  convert  (in) 

21. 

ADC  bit 

0 (out) 

Y. 

E.  0.  C.  status  ADC  (out) 

22. 

Ground 

Z. 

+18  volts  supply 
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Appendix  E contains  the  entire  programming  sequence  for  the  H-8 
computer.  A flow  chart  of  the  program  Is  presented,  then  the  mnemonic 
code  listing  of  the  program,  then  the  definitions  of  the  various 
mnemonic  code  words.  The  actual  machine  language  program  is  then 
listed.  The  parenthesized  numbers  in  the  program  represent  the 
temperature  settings  for  various  rooms.  (11  Is  the  setting  for 
room  A,  (2)  for  room  B,  (31  for  room  C,  (4)  for  room  D,  and  (51 
for  room  I.  These  settings  would  be  programmed  for  the  desired 
temperatures  In  each  room  for  the  time  period  the  temperature  is 
desired  . The  octal  numbers  corresponding  to  the  analog  voltage 
levels  in  the  F.AI-681  are  listed  as  the  last  item  In  this  Appendix. 
The  temperature  levels  associated  with  voltage  levels  are  cxactlv 
ten  times  the  voltage  output  or  one  hundred  time  the  machine  unit 
output  of  the  analog  computer.  Thus  a voltage  of  5.5  volts(.55 
machine  unltsl  represents  55  degrees  F.  The  setting  of  the 
Individual  room  temperatures  would  he  accompl ished  bv  putting  the 
appropriate  octal  number  into  the  appropriate  slot  in  the  program. 
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INITIALIZE 
RFC  IS  t'KRS 


Ann  CONTROL 

V 1)  t O 

OMMANO  WORD 


T r f*  1 1 


INCRF.MF.NT 

RECISTKRS 


OUTPUT 

COMMAS!'  WORD 


INITIALIZE 
COMMAND  WORD 
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MNEMONIC  CODE  PROCKAM 


040  100  MV  1 A 116 

OUT  277 
MV1  A 005 
OUT  277 
MV  I A 116 
OUT  271 
MV I A 005 
OUT  271 

040  120  MV I A 000 
OUT  270 
IN  270 
CPT  060 
JC  300  040 
IN  2 70 
CPI  144 
JC  244  040 
IN  270 
CPT  255 
JC  210  040 

040  151  MV  l A (1) 

STA  100  041 
MV  I A (21 
STA  101  041 
MVI  A (3) 
STA  102  041 
MV l A (4) 
STA  103  041 
MVI  A (5) 
STA  104  041 

-IMP  350  040 

040  210  MVI  A (11 

STA  100  041 
MVI  A (21 
STA  101  041 
MVI  A (31 
STA  102  041 
MVI  A (41 
STA  103  041 
MV  I A (51 
STA  104  041 

JMP  350  040 


040  244 


APPENDIX  E 
MVI  A (1) 

STA  100  041 
MVI  A (2) 

STA  101  041 
MV  1 A (3) 

STA  102  041 
MVI  A (4) 

STA  103  041 
MVI  A (31 
STA  104  041 

JMP  350  040 

040  300  MV  I A (1) 

STA  100  041 
MVI  A (2) 

STA  101  041 
MVI  A (3) 

STA  102  041 
MVI  A (4) 

STA  103  041 
MVI  A (51 
STA  104  041 

JMP  350  040 

040  350  LXI  SP  110  041 
LX I B 064  041 
LXI  D 071  041 
LXI  H 100  041 

040  364  LDAX  B 

OUT'  270 
IN  270 
CMP  M 

.INC  004  041 
XTHL 
LDAX  D 
ORA  M 

STA  076  041 
XTHL 

041  004  LDAX  B 
CPI  005 
JZ  020  041 
INX  B 
I NX  D 
INX  H 

JMP  364  040 
LDA  076  041 
OUT  276 
MVI  A 000 
STA  076  041 
JMP  120  040 


041  020 
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MNEMONIC  CODE 


definitions  7 


27t> 

(Atl*  M (Compute  memotvl 

(A)  ((H)  (I.)) 

The  content  of  llio  memory  location  u’.  ■ ad- 
dress is  contained  in  llio  II  and  I tep.islers  iv 
subtracted  Irotn  the  t ontent  (il  tin1  at  emmd att" 
The  accumulator  remains  unchanged  The  con 
dit ion  flat's  ate  set  as  ,t  result  of  the  sub'!  Van 
The/  this  is  set  t>»  I if  ( \)  ((H)(1))  lh,  i 't  Hat: 

is  set  to  1 if  (A)  < ((H)  (1  )). 


Cycles: 

States: 

Aditressmy 
F lays 


? 

7 

top.  indirect 
/.S.P.CY.AC 


:i?t. 

tiH  data  (Compare  immediate) 

(A)  (byte  2) 

The  t ontent  of  the  set  oml  byte  of  the  histim  lion 
is  subtracted  front  the  content  of  the  a; 
cumulutor  The  condition  flaps  are  set  by  the 
result  of  the  subtrat  turn  The  Z II  a;:  is  set  to  1 it 
(A)  (byte  2).  The  (1Y  tlait  is  set  to  1 it  (A)  (hvtt 
2). 


p 

t , i — ; r , r , t 

rv  \ ' i 

i t i i i 

1 0 1 

1 

L 

do  to  byte 

1 

Cycles 
States 
Atltli  essmy 
Flays 


? 

7 

immediate 

Z.S.P.O  ,AC 


7 Heathktt  Manual  for  the  Digital  Computer  Model  H-8  Oporat  lon_5?s -.201  V 
(Benton  Harbor,  Michigan:  Heath  Company,  1^77)  pastes  40-^7 . 


t 

I 

r 

v 

7 


**  ~ • 


ATPENDIX  . 


CONDITION 

CCC 

NZ 

not  /.ero  (Z  0) 

000 

Z 

zero  (Z  1) 

001 

NC 

no  carry  (CY 

0) 

010 

C 

carry  (CY  1) 

Oil 

PO 

parity  odd  (P 

0) 

too 

Pl\ 

parity  even  (!’ 

11 

101 

I' 

plus  (S  l1) 

1 10 

M 

minus  (S  1) 

m 

.Ill  t 

JMP  addr  (lump) 

(PC)  « 3)  (byte 

Control  is  transferred  to  the  instruction  wluis 
address  is  specified  m byte  2 and  byte  1!  of  ill 
current  instruction. 


P 7 



0 1 o 1 

0 1 0 1 1 

n 

low  order 

addr 

1 

high  Older 

addr 

• 

. . J 

Cycles 

Stales 

Addressing 

Flags 

3 

10 

immediate 

none 

3 (0-7)2 

(condition  .ulilr  (Comlition  jump) 

It  (CCC). 

(PC),  (byte  3)  (byte  2) 

It  the  spm  i fil'd  condition  is  true,  i ontrol  is  tr.ins 
furred  to  the  instruction  whose  address  is 
specified  in  byte  * .mil  liyte  2 ot  the  current  in- 
struction; otherwise,  control  continues  sequen- 
tially. 


r i • i 

L 

c 1 c 1 

1 1 

c 0 1 1 1 0 

low-order 

addr 

high  -order 

addr 

Cycles. 
St. lies: 
Addr  i wing 
Flags 


3 

10 

immediate 

none 


Addressing  d 

f liUjs  none 


003  (ti.C)  0;»3  (H,L) 

0?3  (D,E)  063  (5  ,r) 

\X  rp  (Ini  r.  ment  register  pair) 

(r!i|  (rl)  . (rh)  (r!)  * 1 

t he  1 untent  of  the  register  pair  rp  is  inn  rinentec 
hv  one  N’O'l  Iv  No  ( audition  flaps  are  nffei  ted 


Cycles  1 
States  5 
Addressing  register 
Rags  none 


n?'j 

1. 1)A  atlslr  (load  Aeeunnilatur  direi  ?! 

(A),  ((hvte  it)  (byte  2)) 

the  ( ontent  of  the  memory  lout  I ion.  whose  ,ui 
tlic.s  is  spin  ifictl  ni  hyte  'l  am!  byte  .1  of  th 
instnii  tion.  is  moved  to  the  aci  nine!  0.0' 


tJSH 


R 


t ' 


J 

1 ] 
4 ■ 


f 


(A)  . ((rp» 

The  i mitet  of the  m('iimr\  loiation,  whose  ad 
dross  is  in  (ho  register  jour  rp,  is  nuivod  to  i o ■ stei 
A NOT!'.  Only  register  pairs  rp  1!  (re  . ■ ots  H 
and  O or  rp  0 (resistors  t)  and  !•')  may  l>o 
spo  itiod 

. . J j 1 — - j 

0 1 0 | R 1 F ( 1 0 1 1 1 0 

Cvcics  :> 

States  / 

Altai i sm itij  i vs)  util  .it 

f Ijijs  none 

(l(i  1 ; II,  (.  ) 01  1 (11,1.) 

(>.•  I (11.1)  001  (S.  1’) 

1 \1  rp,  data  It.  (In  d register  pair  nnino.liato' 

(rh)  « (In to  .*i. 

(rl)  s (Into!1) 

Hite  .»  n!  tho  instrm  turn  is  nun  od  into  the  1 ;h 
order  to;;istor  (rh)  o'  tho  register  pair  rp,  Hi : o' 

tlto  mstnii  turn  is  nun  od  into  tho  low  orde'  ; t ops 
tor  (rl'  ot  tho  register  pair  rp 


! 

0 0 

. " ' * I 

o * o 1 

lasv  oulfi 

data 

tinjti  outei 

data 

. 

C voles 

.< 

States 

to 

Adore"  ’ll) 

tinmi'ili.itt- 

l lidjs 

lOlU? 

l)  (0  ')(. 

MM  r,  data 

(Moi o Immediate) 

0(0  "1. 

(r)  . (Into.') 

The  emit  out  ot  In  to  of  tho  in  stun  t .on  is  ntovoi ; 
to  register  i 


r 

r 

L_ 


p^ir; 

da t a hy  t o 


I P 
_l _ 


C voles 
States 
Atklfinsinii 
f (aits 


7 

immediate 

none 
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V M (OH  memory) 

(A).  (A)  v ((11)  (I.)) 

Tin  content  of  tlx;  mem.oi  v location  whose  ad- 
dress is  contained  in  the  U and  !.  registers  is 
im  his  ve-OK'd  with  the  content  of  the  ac- 
cumulator The  result  is  placed  in  the  ar 
curm.'a’ur  The  CY  and  AC  Hags  are  cleared 


Addrcssiny:  icy.  indirect 

Flays  /'.S.P.CY.At 


(Output) 


The  content  of  the  accumulator  is  placed  on  the 
eight-hit  hi-diroetiona!  data  tins  for  transmission 
to  the  specified  port. 


i n du  t 


por  t 


Cycles:  3 

States:  10 

Atidressiny : direct 

Flays:  none 


002 

STA  addr  (Store  Accumulator  direct) 

((byte  3)  (byte  2))  (A) 

The  content  of  the  accumulator  is  moved  to  tl 
memory  location  whose  address  is  specified  i 
byte  2 and  byte  3 of  the  instruction. 


low  order  addr 


high  order  addr 


Addressing 
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H (Kxc  hnngo  stn<  k top  with  I!  a 

(I  ) •>  ((Si'll 

INI  • #((NI*J  < 1) 

1 he  i onlent  nf  tin;  ! i roister  i.s  cm  h r t .1 

the  content  lit  the  memory  location  u . i 

dress  is  spet  ifieil  h\  tile  « ontnit  til  tv 
1 lie  i oMtellt  of  till'  11  : , I’Kster  is  c\,  ) i.i f • ■ . 

the  Content  of  the  mentors  lm  e Oil  \\  . 

itress  ,s  une  more  then  the  i onter.l  of  nr 


T1” 

1~  - r r 

110 

o 1 c 1 1 ' 1 

Cycles 

5 

Suites 

IS 

Atkirejsiiu 

i oil  unliiL-ci 
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OCTAL  CODE  PROGRAM  LISTING 


120 


040  100 

076 

116 

323 

277 

104 

076 

005 

323 

277 

110 

076 

116 

323 

271 

114 

076 

005 

323 

271 

120 

076 

000 

323 

270 

124 

333 

270 

376 

060 

130 

332 

300 

040 

333 

134 

270 

376 

144 

332 

140 

244 

040 

333 

2 ro 

144 

376 

255 

332 

210 

150 

040 

076 

(1) 

062 

154 

100 

041 

076 

(21 

160 

062 

101 

041 

076 

164 

(3) 

062 

' ' 

041 

170 

076 

(41 

062 

103 

174 

041 

076 

(5) 

062 

040  200 

104 

041 

303 

350 

204 

040 

000 

000 

000 

210 

076 

(11 

062 

100 

214 

041 

076 

(21 

062 

220 

101 

041 

07b 

(3) 

224 

062 

102 

041 

07b 

230 

(41 

062 

103 

041 

234 

076 

(51 

062 

104 

240 

041 

303 

350 

040 

244 

07b 

(11 

062 

100 

250 

041 

076 

(21 

062 

254 

101 

041 

076 

(31 

260 

062 

102 

041 

076 

264 

(4) 

062 

103 

041 

270 

076 

(51 

062 

104 

274 

041 

303 

350 

040 

040  300 

076 

(11 

062 

100 

304 

041 

076 

(21 

062 

310 

101 

041 

076 

(3) 

314 

062 

102 

041 

076 

320 

(41 

062 

103 

041 

324 

076 

(5) 

062 

104 

330 

041 

303 

350 

040 

334 

XXX 

XXX 

XXX 

XXX 

340 

XXX 

XXX 

XXX 

XXX 

344 

XXX 

XXX 

XXX 

XXX 

AITKN01X  K 


3S0 

061 

110 

041 

001 

354 

064 

041 

021 

071 

360 

041 

041 

100 

041 

364 

012 

323 

270 

333 

370 

270 

276 

322 

004 

374 

041 

343 

032 

266 

041 

000 

062 

076 

041 

343 

004 

012 

376 

005 

312 

010 

020 

041 

003 

023 

014 

043 

303 

364 

040 

020 

072 

076 

04 1 

323 

024 

276 

076 

000 

062 

030 

076 

041 

303 

120 

034 

040 

XXX 

XXX 

XXX 

040 

XXX 

XXX 

XXX 

XXX 

044 

XXX 

XXX 

XXX 

XXX 

050 

XXX 

XXX 

XXX 

XXX 

054 

XXX 

XXX 

XXX 

XXX 

060 

XXX 

XXX 

XXX 

XXX 

064 

001 

002 

003 

004 

070 

005 

001 

002 

004 

074 

010 

020 

XXX 

XXX 

041 

100 

(n 

(23 

(33 

(43 

104 

(5) 

XXX 

XXX 

XXX 

1 


121 

I 


► 


1 


■RMCFTTSS 


■■■I 


122 
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ANALOG  VOLTAGE  TO  OCTAL  COWERS  TON  TABLE 


9.96 

377 

8.16 

322 

6.36 

245 

4.56] 

170 

2.761 113 

9.92 

376 

8.12 

321 

6.32 

244 

4.52 

167 

2.72 

112 

9.88 

375 

8.08 

320 

6.28 

243 

4.48! 

166 

2.68 

111 

9.84 

374 

8.04 

317 

6.24 

242 

4.44] 

165 

2.64 

HO 

9.80 

373 

8.00 

316 

6.20 

241 

4.40| 

164 

2.60 

' 107 

9.76 

372 

7.96 

315 

6.16 

240 

4.36 

163 

2.56 

106 

9.72 

1 371 

7.92 

314 

6.12 

237 

4.32 

162 

2.52 

! 105 

9.68 

370 

7.88 

313 

6.08 

236 

4.28 

161 

2.48 

104 

9.64 

367 

7.84 

312 

6.04 

,235 

4.24 

160 

2.44 

! 103 

9.60 

366 

7.80 

311 

6.00 

(234 

4.20! 

157 

2.40 

; 102 

9.56 

365 

7.76 

310 

5.96 

233 

4.16 

156 

2.36 

101 

9.52 

364 

7.72 

307 

5.92 

232 

4 . 12 1 

i 155 

2.32 

100 

9.48 

363 

7.68 

306 

5.88 

231 

4.08 

! 154 

2.28 

077 

9.44 

i 362 

7.64 

1305 

5.84 

230 

4.04 

153 

2.24 

076 

9.40 

361 

7.60 

304 

5.80 

227 

4.001 

1 152 

2.20 

075 

9.36 

360 

7.56 

303 

5.76 

226 

3.96 

1151 

2.16 

074 

9.32 

357 

7.52 

302 

5.72 

225 

3.92' 

! 150 

2.12 

073 

9.28 

356 

7.48 

301 

5.68 

224 

3.88, 

147 

2.08 

072 

9.24 

355 

7.44 

300 

5.64 

223 

3.84j 

146 

2.04 

j 071 

9.20 

354 

7.40 

277 

5.60 

222 

3.80; 

145 

2.00 

j 070 

9.16 

353 

7.36 

276 

5.56 

221 

3.76 

144 

1 .96 

|067 

9.12 

352 

7.32 

275 

5.521 

220 

3.72] 

143 

1.92 

066 

9.08! 

351 

7.28 

274 

5.48! 

217 

3.68! 

! 142 

1.88 

065 

9.04! 

350 

7.24 

273 

5.44! 

216 

3.64] 141 

1.841 

064 

9.00] 

347 

7.20 

2 72 

5.40i 

215 

3.60i 140 

1.80 

063 

8.96 

346 

7.16 

271 

5.36! 

214 

3.56] 137 

1.76 

062 

8.92 

345 

7.12 

270 

5.32] 

213 

3.52! 

136 

1.72 

061 

8.88 

344 

7.08 

267 

5.28 

1212 

3.48: 135 

l .68 

060 

8.84 

343 

7.041 

266 

5.24] 

| 211 

3.44! 

134 

1 .64 

(057 

8.80 

342 

7.00l 

265 

5.20! 

210 

3 .40’ 

133 

1.60 

056 

8.76 

341 

6.96 

264 

5.16j 

207 

3.36. 

132 

l.56j 

1 055 

8.72 

340 

6.92] 

263 

5.12| 

206 

3.32' 131 

I.52' 

054 

8.68 

337 

6.88 

262 

5.08 

205 

3 . 28 1 

130 

1 . 48 ! 

053 

8.64 

336 

6.84] 

261 

5.04 

204 

3.24, 

127 

1 .44 

| 052 

8.60 

335 

6.80 

260 

5.00 

203 

3.20 

126 

1 . 40 1 

| 051 

8.56 

334 

6.76' 

257 

4.96 

202 

3.16|  125 

1.36j 

| 050 

8.52 

333 

6.72 

256 

4.92 

201 

3.12! 

124 

1.32 

! 047 

8.48 

332 

6.68 

255 

4.88 

200 

3 .08; 123 

1 . 28  j 

! 046 

8.44 

331 

6.64 

254 

4.84 

177 

3.04' 

122 

1.24 

045 

8.40 

330 

6.60 

253 

4.80 

176 

3.00! 

121 

1 .20 

044 

8.36 

i 327 

6.56 

252 

4.76 

' 175 

2.96 

120 

1.16 

043 

8.32 

i 326 

6.52 

251 

4.72 

! 174 

2.92! 

117 

1.12 

i 042 

8.28 

325 

6.48 

250 

4 .68  i 

j 173 

2.88; 

116 

1 .08 

041 

8.24 

324 

6.44] 

247 

4.64 

1 172 

2.841  115 

1 .04 

040 

8.20 

323 

6.40* 

246 

4.60! 

1171 

2.80!  114 

1 .00 

I 037 

123 
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Appendix  F contains  the  results  of  test  runs  made  in  order  to 
determine  the  optimum  low  temperature  for  an  unused  room.  The  data 
is  presented  such  that  the  horizontal  axis  is  the  temperature  setting 
of  the  unused  room  and  the  vertical  axis  is  the  total  heat  used 
for  the  time  period  being  considered.  The  bar  graphs  at  the  end 
of  t he  Appendix  represent  the  amount  of  heat  used  when  the  rooms 
indicated  are  not  heated.  Runs  were  made  for  a total  time  period 
of  eight  hours  in  all  cases.  The  first  nine  graphs  are  those  resulting 
when  the  unused  room  is  set  to  the  temperature  on  the  x-axis  for  I 
the  entire  time  period.  The  following  graphs  are  the  result  of  the 
unused  room  being  set  to  the  temperature  on  the  x-axis  for  the  time  : 
period  indicated  on  the  graph,  first  3 hours  out  of  eight,  then 
45  minutes  out  of  eight  hours,  then  23  minutes  out  of  eight  hours. 

As  the  graphs  show,  the  trend  is  always  the  same.  The  total  heat 

v 

used  increases  when  the  temperature  setting  of  the  unused  rooni 
is  raised. 
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This  Appendix  contains  the  use  charts  used  to  write  sample 
programs  for  the  digital  controller.  The  first  use  chart  represents 
the  use  habits  of  a family  which  does  not  use  their  house  in  the  day, 
and  does  not  use  the  bedrooms  until  retiring.  The  second  use  chart 
represents  the  use  habits  of  a family  which  does  use  the  house  during 
the  day.  The  charts  represent  typical  habits  for  people  and  do 
not  really  mean  that  there  is  actually  a family  which  follows  these 
use  charts.  The  results  obtained  from  using  the  charts  i6  meant 
to  be  representative  of  what  could  be  achieved  with  the  controller. 
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USE  CHART  1 


HEATED  ' UNHEATEli  j HEATED  j HEATED 


UNCLASSIFIED 
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ROS  (Contlnua  on  rovorao  mUa  II  nacataary  and  IJantllv  by  block  nutnbmr) 


Automatic  control . 

Environmental  engineering  (Buildings) 
Energy  conservation. 

Buildings. 


20  ABSTRACT  (Contlnuo  on  ravaraa  old*  II  nacaaoary  and  Idantlty  b\  block  number) 


This  project  investigated  the  use  of  microprocessors  and  associated 
hardware  to  control  the  heating  of  a building  in  order  to  save  energy. 

A house  was  simulated  on  the  hybrid  computer  and  controlled  by  a mien 
processor-based  digital  controller.  The  work  includes  the  mathematical  model 
simulation,  flow  charts,  computer  programs  for  controller,  etc.  to  use  for 
different  type  of  houses  and  insulations. 

Thests  showed  that  a significant  amount  of  energy  *as  saved  by  using 
the  controller,  as  much  as  30-^0  o/o,  depending  on  the  use-habits  of  the  hou: 


